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Polarized deep-inelastic scattering data on longitudinally polarized hydrogen and deuterium tar- 
gets have been used to determine double spin asymmetries of cross sections. Inclusive and semi- 
inclusive asymmetries for the production of positive and negative pions from hydrogen were obtained 
in a re-analysis of previously published data. Inclusive and semi-inclusive asymmetries for the pro- 
duction of negative and positive pions and kaons were measured on a polarized deuterium target. 
The separate helicity densities for the up and down quarks and the anti-up, anti-down, and strange 
sea quarks were computed from these asymmetries in a "leading order" QCD analysis. The polar- 
ization of the up-quark is positive and that of the down-quark is negative. All extracted sea quark 
polarizations are consistent with zero, and the light quark sea helicity densities are flavor symmetric 
within the experimental uncertainties. First and second moments of the extracted quark helicity 
densities in the measured range are consistent with fits of inclusive data. 

PACS numbers: 13.60.-r, 13.88.+e, 14.20.Dh, 14.65.-q 



I. INTRODUCTION 



Understanding the internal structure of the nucleon re- 
mains a fundamental chaUenge of contemporary hadron 
physics. From studies of deep- inelastic lepton- nucleon 
scattering(DIS), much has been learned about the quark- 
gluon structure of the nucleon, but a clear picture of the 
origins of its spin has yet to emerge. The pioneering 
experiments to explore the spin structure of the nucleon 
performed at SLAC 0,0 were measurements of inclusive 
spin asymmetries, in which only the scattered lepton is 
observed. Until recently, inclusive measurements have 
provided most of the current knowledge of nucleon spin 
structure. The objective of these studies was to deter- 
mine the fraction of the spin of the nucleon which is car- 
ried by the quarks. The nucleon spin can be decomposed 
conceptually into the angular momentum contributions 
of its constituents according to the equation 



(sf) = - = -AI] + L, + J„ (1) 



where the three terms give the contributions to the nu- 
cleon spin from the quark spins, the quark orbital angular 
momentum, and the total angular momentum of the glu- 
ons, respectively. Early calculations based on relativistic 
quark models [3, 0] suggested AS k, 2/3, while more 
precise experiments on DIS at CERN, performed by the 
European Muon Collaboration (EMC) |E0, led to the 
conclusion that AE w 0.1 — 0.2. 

With these indications of the complexity of the spin 
structure, it was quickly realized that a simple leading 
order (LO) analysis that did not include contributions 
from gluons was incomplete. More recent next-to-leading 
order (NLO) treatments provide a picture more appro- 
priate to our present understanding of QCD. The focus 
has been on the polarized structure function gi{x, Q^) for 



the proton, given by 

+ Cs{x,as{Q^))<^AJ:{x,Q^) (2) 
+ 2nfCg{x,asiQ'^)) (E) Ag{x,Q'^)], 

where (e^) — nJ^E"^;^ef , is the electric charge of the 
quark of flavor q, the operator (g) denotes convolution 
over X, Aqns and AS are respectively the nonsinglet and 
singlet quark helicity distributions, and Ag is the gluon 
helicity distribution. Here x is the usual Bjorken scaling 
variable, is the squared four-momentum transfer, 

and nf is the number of active quark flavors. The coef- 
ficient functions Cns, Cs, and Cg have been computed 
up to next-to- leading order d, in . At NLO they 
as well as their associated parton distributions depend 
on the renormalization and factorization schemes. While 
the physical observables are scheme independent, parton 
distributions will be strongly scheme dependent, but con- 
nected from scheme to scheme by well-defined relation- 
ships. In a recent NLO analysis [T3| of available data for 
gi , the SMC group presented results for the first moment 
of gi , which is given by 

j dxg^{x,Q^) = ^[CNs{Q\as{Q''))AqNs{Q'') 
+Cs(g^a,(g2))ao(Q')], 

(3) 

where the dependent quantities C^s, Cs, and Agjvs 
are the first moments over x. In the Adler-Bardeen 
scheme used by the SMC group the singlet axial charge 
flo is 

ao(Q2) = AS-3^^A5(Q2), (4) 

where AS is the first moment of the singlet quark dis- 
tribution, and Ag{Q^) the gluonic first moment. The 
SMC group finds that the analysis of the evolu- 
tion of the world data base gives AS = O.BSlgQg 
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(stat)t0:03(syst)l°:°^(th) and A5(lGeV^) = 0.99^1^1 
(stat)tQ;22(syst)lj!);45(th). The resulting value of the sin- 
glet axial charge is — 0.23 ± 0.07(stat)±0.19(syst). 
While this result strongly constrains the total quark spin 
contribution to the nucleon spin, the limited information 
it provides on the flavor structure of AE is critically de- 
pendent on the assumptions of SU (3) flavor symmetry in 
the interpretation of hyperon beta-decay which are made 
to constrain Aqpfs- A central issue in the analysis of the 
inclusive data from these experiments is their sensitivity 
to SU (3) symmetry breaking, and the reliability of esti- 
mates of the contributions to the first moments coming 
from the unmeasured low x region. 

With rare exceptions, the experiments listed above 
have studied inclusive polarized DIS where only the scat- 
tered lepton is detected. Their sensitivity is limited to 
the polarization of the combination of quarks and anti- 
quarks because the scattering cross section depends on 
the square of the charge of the target parton. The key to 
further progress is more specific probes of the individual 
contributions of Eq. Q to the proton spin. Determina- 
tion of the polarization of the gluons is clearly of high 
priority, and a more precise measurement will eliminate 
a major current ambiguity in the implications of exist- 
ing inclusive data. A more direct determination of the 
strange quark polarization will avoid the need for the use 
of data from hyperon decay and the assumption of SU (3) 
flavor symmetry. Measurements which are sensitive to 
quark flavors will allow the separation of quark and anti- 
quark polarizations. The HERMES experiment attempts 
to achieve these objectives by emphasizing semi-inclusive 
DIS, in which a tt, K, or p is observed in coincidence with 
the scattered lepton. The added dimension of flavor in 
the final hadron provides a valuable probe of the flavor 
dependence and other features of parton helicity distri- 
butions. With the advanced state of inclusive measure- 
ments and the HERMES data with its added dimension 
in the flavor sector, important issues such as measure- 
ments of moments of matrix elements and their accessi- 
bility to measurement can be revisited. Indeed, the re- 
sults reported here, which address the issue of the flavor 
dependence of quark helicity densities, mark the logical 
next step in unraveling the spin structure of the proton. 

This paper begins with a brief development of the for- 
malism required to describe semi-inclusive DIS. It is fol- 
lowed by a description of the HERMES experiment and 
the analysis procedures for flavor tagging which produce 
a comprehensive set of spin asymmetries and a detailed 
flavor decomposition of the quark helicity densities in the 
nucleon. The formalism and experiment are described in 
sections HTl and iHll Sections IIVI and Ivl detail the analysis 
procedures and the resulting cross section asymmetries. 
The extraction of the helicity distributions is explained 
in section IVll while section IVlII summarizes an alterna- 
tive approach to measuring strange quark distributions. 
Partial first and second moments of the extracted he- 
licity distributions and of their singlet and non-singlet 
combinations in the measured kinematic range are given 



in section IVIIII where they are also compared to other 
existing measurements and to results from global QCD 
fits. The conclusions from these results are discussed in 
section IIXI The formalism used for the QED radiative 
and detector smearing corrections is presented in some 
detail in App.fXland tables with the numerical results of 
the present analysis are given in App. ^ 

II. POLARIZED DIS 
A. Polarized Inclusive DIS Formalism 

The main process studied here is depicted in Fig. ^ 
An incoming positron or electron emits a spacelike vir- 
tual photon, which is absorbed by a quark in the nucleon. 
The nucleon is broken up, and the struck quark and the 
target remnant fragment into hadrons in the final state. 
Only the lepton is detected in inclusive measurements 
while detection of one or more hadrons in the final state 
in semi-inclusive measurements adds important informa- 
tion on the scattering process. Contributions from 
exchange can be neglected at the energy of the present 
experiment. 

The kinematic variables relevant for this process are 
listed in Tab. The formalism for DIS is developed in 
many texts on particle physics [iTllT^IT^ . Here, the for- 
malism for polarized DIS is briefly summarized in order 
to introduce the various measured quantities. 




FIG. 1: Diagram of the deep-inelastic scattering process. The 
incoming lepton emits a virtual photon which is absorbed 
by one of the quarks in the nucleon. In the case depicted, 
the struck quark fragments into a pion in the final state. In 
semi-inclusive processes, the scattered lepton and part of the 
hadronic final state are detected in coincidence. 

The inclusive DIS cross-section can be written as fol- 
lows: 

5^ 
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where L^^ is a tensor that describes the emission of the 
virtual photon by the lepton and other radiative pro- 
cesses; it can be calculated in Quantum Electro Dynamics 
(QED). The tensor W'^'^ describes the absorption of the 
virtual photon by the target; it contains all of the infor- 
mation related to the structure of the target. Symmetry 
considerations and conservation laws determine the form 
of W'' {cf. fH fll), which for a spin-1/2 target and 
pure electromagnetic interaction reads: 



Pu 



P-q 



Fi 

Pu^ 



P-q 

'^'^ 



P-q 



(6) 
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TABLE I: Kinematic variables in deep inelastic scattering 

, /Tpi f\ 4-momenta of the initial and 
k = {E, k), k ={E , k') gj^^j g^j^^g leptons 

Q Polar and azimuthal angle of the 

' scattered lepton 

4-momentum of the initial target 
nucleon 

4-momentum of the virtual 
photon 

Negative squared 4-momentum 
transfer 

Energy of the virtual photon 

2p.g - 2Mv Bjorken scaling variable 

Fractional energy of the virtual 
photon 

— {P + q)^ Squared invariant mass of the 

= jV/^ + 2M V — photon-nucleon system 

4-momentum of a hadron in the 
final state 



P (M, 0) 



q = k — k' 



^ _ lab ^ _ ^, 



lab 

y — P-k E 



P = {Eh,p) 



Pq ^ V 



Xp — 



Fractional energy of the observed 
final state hadron 
Longitudinal momentum fraction 
of the hadron 



In this expression, Fi and F2 are unpolarized structure 
functions, while gi and g2 are polarized structure func- 
tions that contribute to the cross section only if both the 
target and the beam are polarized. The usual Minkowski 
metric is given by g^^, and e^ij/a/S is the totally anti- 
symmetric tensor. The four-vector S is the spin of the 
nucleon, and q and P are defined in Tab. In general, 
the structure functions depend on i/ and Q^. They can 
also be defined in terms of the dimensionless scaling vari- 
ables y, the fractional energy transfer to the nucleon, and 
X, the Bjorken scaling variable. The latter is equal to the 



fraction of the nucleon's light-cone momentum carried by 
the struck quark. 

The structure functions are given in the quark-parton 
model (QPM) by: 

(8) 

where the sum is over quark and antiquark flavors, and 
eq is the charge of the quark (or antiquark) in units of 
the elementary charge e. The functions q~^ (x) {q~ (x)) are 
the number densities of quarks or antiquarks with their 
spins in the same (opposite) direction as the spin of the 
nucleon. The structure function Fi{x) measures the total 
quark number density in the nucleon, whereas gi (x) is the 
helicity difference quark number density. Both densities 
are measured as a function of the momentum fraction 
carried by the quark. The structure functions Fi{x) and 
F2{x) are related by the equation 



1 + T 

2xFi(x) = ' F2(x), 



(9) 



which reduces to the well-known Callan-Gross relation 
[1^ in the Bjorken limit. R{x,Q^) is the ratio of lon- 
gi tudinal to transverse DIS cross sections, and 7 = 
y/ jv^ . The structure function 172(2;) vanishes in the 
quark-parton model since it is related to Q suppressed 
longitudinal-transverse interference, which is absent in 
the simple QPM. 

In typical experiments the polarized cross sections are 
not measured directly. Rather, their asymmetry 



Ax = 



^1/2 - 0'3/2 
fl/2 + 0'3/2 



(10) 



is measured, where tTi/2 is the photo-absorption cross 
section for photons whose spin is antiparallel to the tar- 
get nucleon spin, while (73/2 is the corresponding cross 
section for photons whose spin is parallel to the target 
nucleon spin. Angular momentum conservation requires 
that in an infinite momentum frame, the spin-1 photon 
be absorbed by only quarks whose spin is oriented in the 
opposite direction of the photon spin. Consequently, a 
measurement of the difference of these two cross sections 
is related to the polarized structure function 171: 



ffi « a'i/2 - 0-3/2 



(11) 



The structure function F\ is proportional to the sum of 
the cross sections, (T1/2 + cr3/2, with the result that the 
spin structure function g\ can be deduced from A\ by 
using a parameterization of F\ based on world data. 

The picture of the structure functions described to this 
point is based on the quark-parton model of point-like 
constituents in the nucleon. The model can be extended 
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to a more general picture that includes quark interac- 
tions through gluons in the framework of quantum chro- 
modynamics (QCD). In this QCD inspired parton model, 
scaling is violated and the quark densities become de- 
pendent. However, in leading order of the strong coupling 
constant, Eqs. Q and ((HJ still hold if the replacements 
F\(x) — > Fi{x,Q'^) etc. are made. To this order the 
structure functions describe the nucleon structure in any 
hard interaction involving nucleons; they are universal. 



when (72 = 0. This approximation is justified in view of 
the small measur ed values of 32(2;) El El 113 

and the 

kinematic suppression of its contributions in the present 
experiment. The residual effect of the small non-zero 
value of g2{x) is included in the systematic uncertainty 
on Ai as described in section Ivl 



C. Polarized Semi— Inclusive DIS Formalism 



B. Relation to the Inclusive Asymmetries 

While the spin orientation of the nucleon and the vir- 
tual photon is the configuration of primary interest, in 
any experiment only the polarizations of the target and 
the beam can be controlled and measured directly. The 
measured asymmetry Amcas of count rates in the anti- 
aligned and aligned configuration of beam and target po- 
larizations is related to the asymmetry A|| of cross sec- 
tions via 



- Pb Pt fa A 



(12) 



where the kinematic dependencies on x and have been 
dropped for clarity. The factors pb and pt are the beam 
and target polarizations, and fz) is the target dilution 
factor. This quantity fo is the cross section fraction 
that is due to polarizable nucleons in the target (1 for H, 
0.925 for D, and ~ 1/3 for ^He in gas targets; generally 
smaller for other commonly used polarized targets). In 
this experiment the dilution factor fjj is not further re- 
duced by extraneous unpolarized materials in the target 
such as windows, etc. 

The asymmetry in the lepton-nucleon system, A|| , is 
related to the physically significant asymmetry Ai for 
photo-absorption on the nucleon level by 



Au =I?(l + ?77) Ai 



(13) 



where 77 = ejy/[l — (1 — y) e] is a kinematic factor, and 
f/2 ~ is assumed. The factor D depends on x and 
Q^, and accounts for the degree of polarization transfer 
from the lepton to the virtual photon. It is called the 
depolarization factor and is given by 



D = 



1 -l-ei? 



(14) 



where e is the polarization parameter of the virtual pho- 
ton, 



2q\ 20 



1-y 



(15) 



The photon-nucleon asymmetry Ai is related to the 
structure function gi by 



(16) 



As noted in section inclusive polarized DIS is sen- 
sitive only to the sum of the quark and antiquark dis- 
tribution functions because the scattering cross section 
depends on the squared charge of the (anti-)quarks. The 
polarizations of the individual fiavors and anti-flavors are 
accessible in fits to only the inclusive data, where addi- 
tional assumptions are used; e.g. the Bjorken sum rule 
is imposed and the quark sea is assumed to be SU (3) 
symmetric Iil8i] . 

The contributions from the various quarks and anti- 
quarks can be separated more directly if hadrons in the 
final state are detected in coincidence with the scattered 
lepton. Measured fragmentation functions reveal a sta- 
tistical correlation between the flavor of the struck quark 
and the hadron type formed in the fragmentation process. 
This reflects the enhanced probability that the hadron 
will contain the flavor of the struck quark. For example, 
the presence of a 7r+ in the flnal state indicates that it 
is likely that a u-quark or a d-quark was struck in the 
scattering because the 7r+ is a (ud) bound state. The 
technique of detecting hadrons in the final state to iso- 
late contributions to the nucleon spin by specific quark 
and antiquark flavors is called flavor tagging. Note that 
in this case scattering from a u-quark is favored both by 
its charge (2e/3) and by the fact that the d-quark is a 
sea quark and hence has a reduced probability of exist- 
ing in the proton in the x range covered in the analysis 
presented here (0.023 < x < 0.6). 

While the cross section asymmetry Ai is of interest for 
inclusive polarized DIS, the relevant quantity for polar- 
ized semi-inclusive DIS (SIDIS) is the asymmetry in the 
cross sections of produced hadrons in the flnal state: 



'^1/2 



'3/2 



(17) 



'3/2 



in analogy to Eq. H10|l. but where cr'* now refers to the 
semi-inclusive cross section of produced hadrons of type 
h instead of the inclusive cross-section. 

In analogy to Eq. Q , the semi-inclusive DIS cross 
section can be written as 



dV 



dx dQ^ dz dp^ 1 



oc 



(18) 



where the hadron tensor, Wj^" now contains additional 
degrees of freedom corresponding to the fractional en- 
ergy z of the flnal state hadron, the component px of the 
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final hadron three momentum transverse to that of the 
virtual photon, and the azimuthal angle 4> of the hadron 
production plane relative to the lepton scattering plane. 
Integration over (j) and produces the cross section rel- 
evant to the present experiment. The assumption of fac- 
torization permits the separation of the hadron degrees 
of freedom from the variables associated with the lepton 
vertex. Consequently, kinematic factors depending only 
on X and Q^, e.g. the depolarization factor D, are car- 
ried over directly from inclusive scattering in relating the 
semi- inclusive asymmetry Ay to A\. 

In leading order, the resulting cross sections <J^f2 ^^'^ 
^3/2 written in terms of the quark distribution 

functions and fragmentation functions D^{z, Q^): 



/2(3/2) 



dx dQ^ dz 



q+^->{x,Q')D^{z,Q'), (19) 



where the dependence on the kinematics is made explicit. 
The fragmentation function is a measure of the prob- 
ability that a quark of flavor q will fragment into a hadron 
of type h. 

A procedure identical to that described in the previ- 
ous subsection relates the measured quantity ^moas to 
the photon-nucleon asymmetry A^. The latter can be 
expressed in terms of the quark helicity densities Aq and 
the fragmentation functions: 



^^,^^^^,^^^_j:,elAqix,Q^)D^iz,Q-) 



J:,, 4 q'{x,Q^)D>^,{z,Q^) 
This equation can be rewritten as follows: 



Al{x,Q^,z) = J2'Pgix,Q',^) 



q{x,Q^) 



(20) 



(21) 



where the quark polarizations (Aq/q) are factored out 
and purities are introduced. The purity is the con- 
ditional probability that a hadron of type h observed in 
the final state originated from a struck quark of flavor 
q in the case that the beam/target is unpolarized. It is 
related to the fragmentation functions by: 



eU{x,Q^)D^{z,Q^) 



E,' e^, q'ix,Q^ 



^D'Mz,Q^) 



(22) 



This concept of a purity is generalized to inclusive scat- 
tering by setting the fragmentation functions to unity in 
Eq. H22|) . This allows the inclusion of the inclusive data 
in the same formalism. 

The determination of the quark polarizations using fla- 
vor tagging based on Eq. trades the assumptions 
used in global fits to inclusive data for the modeling of 
the fragmentation process. The purity formalism based 
on Eq. H19|) additionally implies the factorization of the 
hard scattering reaction and the fragmentation process. 



In the analysis presented here the purities were calcu- 
lated from a Monte Carlo simulation of the entire scat- 
tering process. The determination of the purities and 
the extraction of the quark polarizations on the basis of 
Eq. 121(1 are explained in detail in section IVTI 



III. EXPERIMENT 

The HERMES experiment is located in the East Hall 
of the HERA facility at DESY (see Fig. 121 . Although 
HERA accelerates both electrons (or positrons) and pro- 
tons, only the lepton beam is used by HERMES in 
a fixed-target configuration. The proton beam passes 
through the mid-plane of the experiment. The target is 
a gas cell internal to the lepton ring. There are three 
major components to the HERMES experiment: the po- 
larized beam, the polarized target, and the spectrometer. 
All three are described in detail elsewhere. As this paper 
reports on data collected from the years 1996 until 2000, 
the following describes the experimental status during 
this time. 



A. Polarized Beam 

Detailed descriptions of the polarized beam, the beam 
polarimeters, and the spin rotators are given in Refs. [l9l 
lis m 123. The electron/positron beam at HERA is 
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FIG. 2: Schematic diagram of the HERA accelerator layout 
until 2000 with the location of the four experiments. Also 
shown are the locations of the spin-rotators and the two po- 
larimeters. 

self-polarized by the Sokolov-Ternov mechanism [23l |. 
which exploits a slight asymmetry in the emission of 
synchrotron radiation, depending on whether the spin 
of the electron/positron in the spin flip associated with 
the emission is parallel or anti-parallel to the magnetic 
gui de field. This very small asymmetry (one part in 10^*^ 
'24'|), causes the polarization of the beam to grow asymp- 
totically, with a time constant that depends on the final 
polarization of the beam. This time constant is used in 
dedicated runs to verify the calibration of the polarime- 
ters that measure the degree of polarization of the beam. 

Polarizing the lepton beam at HERA is therefore a 
matter of minimizing depolarizing effects rather than one 
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of producing an a priori polarized beam and keeping it 
polarized. An unpolarized beam is injected into the stor- 
age ring and polarization builds up over time, typically 
in 30-40 minutes. 

The Sokolov-Ternov mechanism polarizes the beam in 
the transverse direction, i.e. the beam spin orientation is 
perpendicular to the momentum. The beam spin orien- 
tation is rotated into the longitudinal direction just up- 
stream of HERMES, and is rotated back into the trans- 
verse direction downstream of the spectrometer. The lo- 
cations of the spin rotators are indicated in Fig. 13 

The beam polarization is measured continuously by 
two instruments, both based on asymmetries in the 
Compton backscattering of polarized laser light from the 
lepton beam. The transverse polarimeter |20l |2]| , mea- 
sures the polarization of the lepton beam at a point where 
it is polarized in the transverse direction. The interac- 
tion point (IP) of the polarized light with the lepton beam 
is located about 120 m downstream of the HERA West 
Hall (see Fig. The polarimeter uses a spatial (up- 
down) asymmetry in the back-scattering of laser light 
from the polarized lepton beam. Back-scattered pho- 
tons are measured in a split lead-scintillator sampling 
calorimeter, where the change in the position of the pho- 
tons with initial circular polarization determines the po- 
larization of the lepton beam. The calorimeter is located 
70 m downstream of the IP. 

A second polarimeter some 90 m downstream of 
the HERMES target (and just before the spin is rotated 
back to the transverse direction) measures the polariza- 
tion of the beam when it is in the longitudinal orienta- 
tion. This polarimeter is also based on Compton back- 
scattering of laser light, but in this case the asymme- 
try is in the total cross-section, and not in the spatial 
distribution. The larger asymmetry in this case allows 
a more precise measurement of the beam polarization. 
The higher precision is reflected in the smaller system- 
atic uncertainties of the polarization measurements in the 
years 1999 and 2000. Additionally, this second polarime- 
ter provides the possibility to measure the polarization 
of each individual positron bunch in HERA. This feature 
is particularly useful for the optimization of the beam 
polarization when the HERA lepton beam is in collision 
with the HERA proton beam. The existence of two po- 
larimeters also allows a cross-check of the polarization 
measurement to be made. 

The beam polarization was typically greater than 50 % 
in the later years of the experiment, attaining values near 
60 % for many fills of the storage ring. Average beam 
polarizations, the precision of the polarization measure- 
ment, as well as the charge of the HERA lepton beam for 
each year are given in Tab. ^ 



B. Polarized Target 

HERMES has used two types of polarized targets over 
the years. In 1995 an optically pumped polarized '^He 



TABLE IL Beam polarization and HERA lepton beam charge 
for each year of HERMES running covered in this paper. The 
numbers are weighted by the luminosity so that the value 
at the beginning of the fill dominates. Polarization values 
were larger at the end of the fill. The fractional uncertainties 
quoted are the ones used in the data analysis. 



Year 


Lepton 
beam charge 


Average 
Polarization 


Fractional 
Uncertainty 


1996 


e+ 


52.8 % 


3.4% 


1997 


e+ 


53.1% 


3.4% 


1998 


e~ 


52.1% 


3.4% 


1999 


e+ 


53.3 % 


1.8% 


2000 


e+ 


53.3 % 


1.9% 



target was installed • Since these data are not used in 
the present analysis, no description of this target is given 
here. In 1996-97, polarized hydrogen was used, while in 
1998-2000 the target was polarized deuterium. In both 
cases, the source of polarized atoms was an atomic beam 
source (ABS). The ABS and the Breit-Rabi polarime- 
ter (BRP) used to monitor the degree of polarization are 
described in [2^ |23, H^. A schematic diagram of the 
polarized target is shown in Fig. |3| Briefly, the atomic 
beam source is based on the Stern-Gerlach effect. Neu- 
tral atomic hydrogen or deuterium is produced in a dis- 
sociator and is formed into a beam using a cooled nozzle, 
collimators and a series of differential pumping stations. 
A succession of magnetic sextupoles and radio-frequency 
(RF) fields are used to select one (or two) particular 
atomic hyperfine states that have a given nuclear polar- 
ization. 

The ABS feeds a storage cell which serves to in- 
crease the density by two orders of magnitude. This 
storage cell is located in the HERA ring vacuum and 
it is cooled to a temperature between 70 K (deuterium) 
and 100 K (hydrogen). The cell is 40cm long and had 
elliptical cross-sectional dimensions of 29 x 9.8 mm^ in 
1996-1999 and 21 x 8.9 mm^ in 2000. 

The polarization and the atomic fraction of the target 
were monitored by sampling the gas in the target cell 
using the target gas analyzer (TGA) and the BRP. The 
TGA is a quadrupole mass spectrometer (QMS) which 
measures the relative fluxes of atomic and molecular hy- 
drogen or deuterium and thereby determines the molecu- 
lar fraction of the target gas. The BRP works essentially 
in reverse to the ABS. Single atomic hyperfine states are 
isolated using magnetic and RF fields and the atoms in 
each hyperfine state are counted using again a QMS. The 
electromagnetic fields are varied in a sequence such that 
atoms in each hyperfine state are counted in succession. 
More details are given in [sS |3l| . Tab. Hn lists 

the target type, average polarization, and uncertainty for 
each data set. The differences in the systematic uncer- 
tainties are largely due to varying running conditions and 
the quality of the target cell in use. 
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FIG. 3: Diagram of the HERMES polarized target. Shown are the atomic beam source (ABS), the target gas analyzer (TGA) 
and the Breit-Rabi polarimeter (BRP). SET, MET, and WET label the strong, medium, and weak field transitions in the ABS 
and the BRP. 



TABLE HI: Target type and polarization for each year of 
HERMES running. These values are weighted by the lumi- 
nosity. The uncertainty quoted is the one used in the data 
analysis. 







Average 


Eractional 


Year 


Typo 


Polarization 


Uncertainty 


1996 


H 


75.9 % 


5.5% 


1997 


H 


85.0 % 


3.8% 


1998 


D 


85.6 % 


7.5% 


1999 


D 


83.2 % 


7.0% 


2000 


D 


+85.1,-84.0% 


3.5% 



C. The HERMES Spectrometer 

The HERMES spectrometer is described in detail in 
[3^ . It is a forward spectrometer with large acceptance 
that can detect the scattered electron/positron as well 
as hadrons in coincidence. This allows semi-inclusive 
measurements of the polarized DIS process, which are 
the focus of this paper. A diagram of the spectrometer 
is shown in Fig. 0] 

Briefly, the HERMES spectrometer consists of multi- 
ple tracking stages before and after a 1.3 Tm dipole mag- 
net, as well as extensive particle identification. The ge- 
ometrical acceptance of the spectrometer is ±170mrad 
in the horizontal direction and between ±40 mrad and 
±140mrad in the vertical. The range of scattering an- 
gles is therefore 40 mrad to 220 mrad. The spectrometer 
is split into two halves (top/bottom) due to the need for 
a flux shielding plate in the mid-plane of the magnet 
to eliminate deflection of the primary lepton and proton 
beams, which pass through the spectrometer. The lepton 
beam passes along the central axis of the spectrometer. 



The proton beam traverses the spectrometer parallel to 
the lepton beam but displaced horizontally by 71.4 cm. 



1. Particle Tracking 

The tracking system serves several functions: 

• Determine the event vertex to ensure the event 
came from the target gas, not from the walls of 
the target cell or from the collimators upstream of 
the target. 

• Measure the scattering angles of all particles. 

• Measure the particle momentum from the deflec- 
tion of the track in the spectrometer magnet. 

• Identify hits in the PID detectors associated with 
each track. 

The tracking system consists of 51 planes of wire cham- 
bers and six planes of microstrip gas detectors. Because 
of the width of the tracking detectors in the rear section 
of the spectrometer, it was not possible to use horizon- 
tal wires in these chambers. Instead, wires tilted ±30° 
from the vertical were used {U and V planes), together 
with vertical wires {X planes). All chambers have this 
geometry to simplify the tracking algorithm (see below). 

The majority of the tracking detectors are horizontal 
drift chambers with alternating anode and cathode wires 
between two cathode foils. The chambers are assem- 
bled in modules of six layers in three coordinate doublets 
{XX', UU', and VV). The primed planes are offset by 
a half-cell to resolve left -right ambiguities. 
In order starting at the target, the tracking chambers 
are: 
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FIG. 4: Side view of the HERMES spectrometer. The positron beam enters from the left. The spectrometer is spht into two 
halves, one above the beam and one below, by a flux exclusion plate to protect the beams from the magnetic field. See the text 
for further details on the detectors. 



a. Vertex Chambers (VCl/2): The main purpose of 
the vertex chambers '33l l34| is to measure the scattering 
angle to high precision and determine the vertex posi- 
tion of the interaction. Because of severe geometrical 
constraints and the high flux of particles in the region 
so close to the target, microstrip gas chambers (MSGC) 
were chosen for the VCs. Each of the upper and lower VC 
detectors consist of six planes grouped into two modules 
(VUX and XUV for VCl and VC2 respectively). The 
pitch of the strips is 0.193 mm. 

6. Drift Vertex Chambers (DVC): The drift vertex 
chambers have a cell size of 6 mm in each module in the 
XX'UU'W geometry. 

c. Front Chambers (FCl/2): The front chambers 
[35| are drift chambers mounted on the front face of 
the spectrometer magnet. The cell size is 7 mm in the 
XX'UU'W geometry. 

d. Magnet Chambers (MCI -3): The magnet cham- 
bers j36l| are located in the magnet gap. The MCs are 
proportional wire chambers with a cell width of 2 mm. 
Each module consists of three planes in the UXV geom- 
etry. 

e. Back Chambers (BCl-4)- The back chambers 
[3^ are large drift chambers located behind the spec- 



trometer magnet. The cell width is 15 mm in the 
UU'XX'W geometry. 

A tracking algorithm [s^ defines tracks in front of and 
behind the magnet and the momentum of the scattered 
particles can therefore be determined. The MSGCs con- 
tained in the VCs were not available after 1998 due to 
radiation damage. In their place, vertex determination 
was accomplished by a refined tracking algorithm that 
used data from the FCs together with the point defined 
by the the intersection of the track in the rear of the 
spectrometer (the back-track) with the mid-plane of the 
magnet as an additional tracking parameter. The track- 
ing algorithm is described in more detail in section HVI 

Note that the magnet chambers are used only to track 
particles that do not reach the back of the spectrometer. 
They are useful for the measurement of partial tracks 
(mostly low-energy pions) that can, under certain con- 
ditions, increase the acceptance for the reconstruction of 
short-lived particles, such as A particles. However, these 
chambers as well as the vertex and the drift vertex cham- 
bers are not used in the analysis reported in this paper. 

Multiple scattering, and bremsstrahlung in the case of 
electrons or positrons, in the windows and other detector 
and target cell material which the particle tracks traverse 
limit the momentum resolution of the spectrometer. Af- 
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ter its installation in 1998, the RICH detector because 
of its aerogel radiator assembly and heavy gas radiator 
increased this limit significantly. Plots of the momentum 
and angular resolution are shown in section Hvl 



Z. Particle Identification 

There are several particle identification (PID) detec- 
tors in the HERMES spectrometer. Electrons and po- 
sitrons are identified by the combination of a lead-glass 
calorimeter, a scintillator hodoscope preceded by two ra- 
diation lengths of lead (the pre-shower detector), and 
a transition-radiation detector (TRD). A Cerenkov de- 
tector was used primarily for pion identification. The 
threshold detector was replaced by a Ring-Imaging 
Cerenkov (RICH) detector in 1998. The RICH al- 
lowed pions, kaons, and protons to be separated. Both 
Cerenkov detectors also helped in lepton identification. 

a. The Calorimeter: The calorimeter [39^ has the 
following functions: suppress hadrons by a factor of 10 
in the trigger and 100 offline; measure the energy of elec- 
trons/positrons and also of photons from other sources, 
e.g. 7r° and r\ decays. It consists of two halves each con- 
taining 420 blocks (42 x 10) of radiation resistant FlOl 
lead-glass. The blocks are 9x9 cm^ by 50 cm deep 
(about 18 radiation lengths). Each block is viewed from 
the back by a photomultiplier tube. 

The response of the calorimeter blocks was studied in 
a test beam with a 3 x 3 array. The response to electrons 
was found to be linear within 1 % over the energy range 
1-30 GeV. The energy res olution w as measured to be 
a(E)lE [%] = (5.1 ± l.l)/V£^[GeV] -h (1.5 ± 0.5) 

h. The Pre-Shower Detector: The calorimeter is 
preceded by a scintillator hodoscope (H2) that has two 
radiation lengths of lead in front of it. The hodoscope H2 
therefore acts as a pre-shower detector and contributes 
to the lepton identification. This detector consists of 42 
vertical scintillator modules in each of two halves. Each 
paddle is 1 cm thick and 9.3 x 91 cm^ in area. 

The lead preceding the hodoscope initiates showers 
for leptons but with a much reduced probability for 
hadrons. Pions deposit only about 2 MeV of energy 
on average while electrons/positrons deposit roughly 20- 
40 MeV. H2 suppresses hadrons by a factor of about 10 
with 95 % efficiency for detection of electrons/positrons. 

c. Transition Radiation Detector: The transition 
radiation detector (TRD) rejects hadrons by a factor 
exceeding 300 at 90% electron/positron detection 
efficiency. Each of the upper and lower halves of the 
spectrometer contains six TRD modules with an active 
area of 325 x 75 cm^. Each module consists of a radiator 
and a proportional wire chamber to detect the TR 
photons. The radiators consist of a pseudo-random but 
predominantly two-dimensional array of polyethylene 
fibers with 17-20 /im diameter. The proportional cham- 
bers have a wire spacing of 1.27cm, use Xe:CH4 (90:10) 
gas, and are 2.54 cm thick. 



d. Cerenkov Detector: In 1995-97, a threshold 
Cerenkov detector was operated, which was located be- 
tween the two sets of back tracking chambers. During the 
1996 and 1997 data taking periods, a mixture of 70 % ni- 
trogen and 30 % C4F10 was used as the radiator, resulting 
in momentum thresholds for pions, kaons, and protons of 
3.8, 13.6, and 25.8 GeV respectively. 

As for the other components of the spectrometer, the 
Cerenkov detector consists of two identical units in the 
upper and lower half of the spectrometer. The numbers 
given in the following refer to one detector half. An ar- 
ray of 20 spherical mirrors (radius of curvature: 156 cm) 
mounted at the rear of the gas volume focused the 
Cerenkov photons onto phototubes of diameter 12.7 cm. 
Hinterberger-Winston light cones with an entrance di- 
ameter of 21.7 cm helped maximize light collection. The 
mean number of photoelectrons for a /3 « 1 particle was 
measured to be around five. 

e. RICH: The threshold Cerenkov detector was re- 
placed in 1998 by a ring imaging Cerenkov detector 
(RICH) which allowed kaons and protons to be identified 
as well as pions 40] . The RICH uses a novel two-radiator 
design to achieve separation of pions, kaons, and protons 
over the entire kinematic range of interest (4-13.8 GeV; 
see Fig. lion . One of the radiators is C4F10 gas with an in- 
dex of refraction of n = 1.0014, while the second radiator 
consists of aerogel tiles with index of refraction n — 1.03 
mounted just behind the entrance window. The aerogel 
tiles are 1.1cm thick and they are stacked in five layers 
for a total length of 5.5 cm. A mirror array with a ra- 
dius of curvature of 220 cm focuses the Cerenkov photons 
onto 1934 photomultiplier tubes of 1.92 cm diameter per 
detector half. Details on the analysis of the RICH data 
are given in section Hvl 

/. PID Detector Performance: Plots of the re- 
sponses of the PID detectors are shown in Fig. |5| A 
description of the PID analysis, integrating all the 
detectors, is given in section Hvl 



3. Event Trigger 

Before discussing the trigger itself, two more detectors 
used specifically for the trigger must be introduced: two 
scintillator hodoscopes HO and HI. The hodoscope HI is 
identical to H2 except that it does not have any lead in 
front of it. It is situated between BC4 and the TRD. The 
scintillator HO was added after the first year of running 
to help discriminate against particles traveling backwards 
in the spectrometer. These particles originate in show- 
ers initiated by the proton beam. The HO hodoscope 
is placed just in front of the magnet and therefore has 
enough separation from HI and H2 that it can determine 
whether a particle is going forward or backwards in the 
spectrometer. 

The DIS trigger selects electron/positron events by re- 
quiring hits in the three scintillator hodoscopes (HO, HI, 
and H2) together with sufficient energy deposited in two 
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FIG. 5: Typical PID detector responses. The distributions are based on a small set of the data collected in 2000, except for 
the threshold Cerenkov response, which was computed from a data set of similar size collected in 1997. The truncated mean 
is shown in the case of the TRD. The relative size of the lepton (dashed line) and hadron distributions (solid line) was scaled 
to the flux ratio in the respective data-taking periods to give a better idea of the level of contamination possible from each 
detector. The flux ratio of electrons to hadrons is typically ~ 10 % for these data. 



adjacent columns of the calorimeter, in coincidence with 
the accelerator bunch signal (HERA clock) . The require- 
ment of hits in HO and HI suppresses neutral particle 
background. The calorimeter has a high efficiency for 
electromagnetic showers, but relatively low efficiency for 
hadronic showers. The calorimeter threshold was set at 
1.4 GeV (3.5 GeV for the first period in 1996). 



4- Luminosity Monitor 

The luminosity was measured using elastic scattering 
of beam particles by the electrons in the target gas: 
Bhabha scattering and annihilation for a positron beam, 
M0ller scattering for an electron beam [IJ. The scat- 
tered particles exit the beam pipe 7.2 m downstream of 
the target. They are detected in coincidence by a pair of 



small calorimeters with a horizontal acceptance of 4.6- 
8.9 mrad. The calorimeters consist of Cerenkov crystals 
of NaBi(W04)2 that are highly resistant to radiation 
damage. 



5. Data Acquisition System and Event Structure 

The backbone of the data acquisition system is con- 
structed in Fastbus. It consists of 10 front-end crates, 
an event collector crate, and an event receiver crate, con- 
nected to the online workstation cluster via two SCSI 
interfaces. CERN Host Interfaces (CHI) act as Fastbus 
masters, and their performance is enhanced by Struck 
Fastbus Readout Engines (FRE) containing two Mo- 
torola 96002 DSPs. 

The drift chambers were read out by LeCroy multi- 
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hit, multi-event 16-bit 96 channel TDCs (model 1877). 
Charge from the photomultipliers and from the TRD was 
digitized by LeCroy multi-event 64 channel 1881M multi 
block ADCs. These ADCs and the TDCs are capable of 
sparsifying the data, i.e. online suppressing channels with 
pedestal levels from the readout. The magnet chamber 
readout was instrumented with the LeCroy VME based 
PC0S4 system. The vertex chamber data arrived from 
the detector as a 16 bit ECL STR330/ECL data stream 
and were processed in one of the VC DSPs. Double 
buffering was implemented in the dual DSPs of the Fast- 
bus masters. Event collection on one DSP was done in 
conjunction with readout from the second DSP to the 
DAQ computer. 

In addition to the standard readout, a series of asyn- 
chronous independent events from the luminosity moni- 
tor and from monitoring equipment could be read out at 
rates exceeding 5 kHz. One VME branch with 4 crates 
and three CAMAC branches with 9 crates were used for 
these events. The DAQ dead-time was typically less than 
10% with a total trigger rate of about 300 Hz. 

The data are arranged into the following time struc- 
ture: 

• Burst: Events are grouped into bursts, defined as 
the interval between two successive reads of the ex- 
periment scalers. A burst is roughly 10 s long. Data 
quality is checked on the burst level. 

• Run: The size of the files stored on disk and tape is 
adjusted so that an integral number of runs can fit 
on a tape. At high instantaneous luminosity, one 
run can be as short as lOmin. A run is the basic 
unit of data for analysis. Calibration constants are 
applied at the run level, although not all detectors 
are calibrated with this time granularity. 

• Fill: Runs are grouped into fills, which are simply 
defined as data collected during a given fill of the 
electron storage ring. 

IV. DATA ANALYSIS 

A. Data Quality 

The data used to compute the asymmetries and mul- 
tiplicities were selected by a number of quality criteria 
applied at the burst level: 

• the beam polarization was between 30 % and 80 %, 
and the beam current was between 5 and 50 mA. 
(The upper bounds are beyond values observed dur- 
ing data-taking. They are imposed to reject faulty 
records.), 

• the trigger dead-time was less than 50% and the 
data acquisition worked satisfactorily, 

• the PID system and the tracking detectors worked 
properly. 



• there were no high voltage trips in any of the de- 
tectors, 



• the experiment was in polarized running mode, 

and the target system was required to be fully opera- 
tional. This requirement resulted in polarizations in ex- 
cess of 75 % for both the hydrogen and deuterium targets. 
See Refs. 113,53 more details on the selection of good 
quality data. 



B. Tracking Algorithm 

Particle tracks were reconstructed using the pattern of 
hits in the front and back tracking systems [Sg- In the 
first step of this procedure, the partial front and back 
tracks, which are approximately straight lines, are re- 
constructed separately in each of the U, V, and X ori- 
entations. The algorithm is based on a fast tree search. 
For each orientation, the algorithm begins by considering 
the entire plane and successively doubles the resolution 
by discarding the halves without a hit. In each step the 
combined patterns of all planes in a given orientation 
are compared to a data base of physically possible tracks 
and only corresponding patterns are kept. After about 
11 steps the search reaches a resolution that is sufficient 
for track finding. The projections in the three planes are 
then combined to form the partial tracks in the front and 
the back respectively. 

The front and back tracks are associated by matching 
pairs that intersect in the center of the magnet within 
a given tolerance. For each associated pair, the front 
track is forced to agree with the magnet mid-point of 
the back track, and the front track is recomputed ac- 
cordingly. This procedure improves the resolution of the 
front tracking system, which relies on the EC chambers 
since only they were installed and operational during the 
entire data taking period from 1996 until 2000. Because 
the tracking information from the other chambers was 
not available for this entire period, they were not used 
in order to avoid possible biases for different data tak- 
ing periods. The particle momentum is determined us- 
ing another data base of 520, 000 tracks which contains 
the momentum as a function of the front and back track 
parameters. Multiple scattering in the spectrometer ma- 
terial leads to reduced resolutions of less than 0.03 for the 
reconstructed track momenta and less than 1.5mrad for 
the reconstructed scattering angles. Fig.lHlshows the res- 
olutions for the deuterium data sample as obtained from 
a Monte Carlo simulation of the entire spectrometer. The 
momentum and angular resolution of the hydrogen data 
are better, because of the shorter radiation length of the 
Cerenkov detector compared to the RICH. 
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FIG. 6: Tracking system resolution for lepton and hadron 
tracks for the detector configuration used since 1998. In the 
left panel the relative momentum resolution is displayed, and 
the right panel shows the resolution in the horizontal scat- 
tering angle 6x, both as a function of the track momentum 
P- 



C. Particle Identification Algorithm 

The PID system discriminates between electrons/po- 
sitrons (referred to as leptons in the following), pions, 
kaons, and other hadrons. It provides a factor of about 
10 in hadron suppression at the trigger level to keep data 
acquisition rates reasonable. The hadron rate from photo 
production exceeds the DIS rate by a factor of up to 
400:1 in some kinematic regions. In offline analysis, the 
HERMES PID system suppresses hadrons misidentified 
as leptons by as much as 10* with respect to the total 
number of hadrons, while identifying leptons with effi- 
ciencies exceeding 98%. 

The identification of hadrons and leptons is based on 
a Bayesian algorithm that uses the conditional proba- 
bility P{A\B) defined as the probability that A is true, 
given that B was observed. For each track the condi- 
tional probability P{Hii^i^)\E,p, 0) that the track is a lep- 
ton (hadron) is calculated as 



j:^=iMP{mP,0) P{E\H,,p) 



(23) 



Here is the hypothesis that the track is a lepton 

(hadron), E the response of the considered detector, and 
p and 9 are the track's momentum and polar angle. The 
parent distributions P{E\Hi(^}^),p) of each detector (i.e. 
the typical detector responses) were extracted from data 
with stringent restrictions on the other PID detectors to 
isolate a particular particle type. See Fig. [S] for plots of 
the individual PID detector responses. 

In a first approximation, uniform fluxes P{Hi\p,9) = 
P{Hfi\p, 9) are assumed so that the ratio 



log 
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P{Hi\E,p,t 
P{Hh\E,p,t 



reduces to: 



PID, 



det 



logi 



P{E\Hup) 
PiE\Hh,p) 



(24) 



(25) 



The quantity PID^et is defined for the calorimeter 
(caZ), the pre-shower detector (pre), the Cerenkov detec- 
tor (cer) (the RICH detector (ric) since 1998), and the 
TRD (trd). In the case of the RICH and the TRD this 
ratio is the sum over the PID values of the two radiators 
and the six TRD modules respectively. The PID distri- 
bution of the TRD (PID5) is shown in Fig. □ versus the 
sum of the PID values of the calorimeter, the pre-shower, 
and the threshold Cerenkov/RICH (PID3). The leptons 
(small bump) are seen to be clearly separable from the 
hadrons (large peak). 



xlO' 

0) 

c 3500 
O 3000 
2500 
2000 
1500 
1000 
500 




FIG. 7; Two-dimensional distribution of PID values for all 
particles in the acceptance. The quantities PID3 and PID5 
are defined as PID3 = PIDca! + PIDp^e + PID^ic and PID5 = 

Pl'Dtrd- 



The particle fluxes P{Hni^-^\p,9) were computed in an 
iterative procedure by comparing the calculated ratio 
Eq. H24|l to data and varying the fluxes. These fluxes 
were then combined with PID3 -f PID5 to form the total 
PID value 



PID = PID3 + PID5 - logic * 



(26) 



where $ = P(Hii\p,9)/P{Hi\p,9) is the ratio of hadron 
and lepton fluxes. A plot of the quantity PID which 
was used to discriminate hadrons and leptons is shown 
in Fig. IHl where the two peaks for hadrons and leptons 
are seen to be well separated. Hadrons and leptons were 
identified with limits requiring PID < and PID > 1 re- 
spectively. The lepton restriction provided excellent dis- 
crimination of DIS leptons from the large hadron back- 
ground, with efficiencies larger than 98 % and contamina- 
tions below 1.0% over the entire range in x (see Fig.^. 
Semi-inclusive hadrons were identified with efficiencies 
larger than 99 % and lepton contaminations smaller than 
1.0%, which were determined 42] from data collected 
during the normal operation of the experiment. 
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FIG. 8: The distribution of the total PID value. This loga- 
rithmic ratio of probabilities includes the particle fluxes and 
the responses of all PID detectors. The left hand peak is the 
hadron peak, while the right hand peak originates from lep- 
tons. The limits that were applied in the analysis are shown 
as vertical lines. 
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FIG. 9: Identification efficiency and hadron contamination of 
the DIS lepton sample as a function of x. Because correlations 
between the responses of the PID detectors were neglected, 
the contaminations are uncertain by a factor of two. The 
deuteron data have slightly worse efficiencies and contamina- 
tions because of the better hadron-lepton discrimination of 
the threshold Cerenkov counter compared to the RICH. 



tency with the threshold Cerenkov detector. The pattern 
of Cerenkov photons emitted by tracks passing through 
the aerogel or the gas radiators on the photomultipher 
matrix was associated with tracks using inverse ray trac- 
ing. For each particle track, each hadron hypothesis, and 
each hypothesis for the radiator emitting the photons, 
aerogel or gas, the photon emission angle was computed. 
The average Cerenkov angles {0)°^'% ^ were calculated for 
each radiator {a,g) and particle hypothesis {tt,K,p) by 
including only photons with emission angles within 2 ag 
about the theoretically expected emission angle d^^^xp'^ i 
where ere ~ 8 mrad is the single photon resolution. This 
procedure rejects background photons, and photons due 
to other tracks or the other radiator. Fig. 1101 shows the 
distribution of angles in the two radiators as a function 
of the particle momentum. Based on the Gaussian like- 
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FIG. 10: Cerenkov angles associated with the three parti- 
cle hypotheses as a function of the particle momentum. The 
characteristic angles of Cerenkov light emitted in the aerogel 
(n = 1.03) are given by the solid lines. The characteristic 
angles for emission in the gas (n — 1.0014) are shown as the 
dashed lines. The corresponding histogram entries are exper- 
imentally determined angles of a sample of SIDIS hadrons. 

lihood, 



D. The Cerenkov detectors and hadron 
identification 



The threshold Cerenkov detector identified pions with 
momenta between 4 and 13.8 GeV. A hadron track was 
identified as a pion, if the number of detected photo- 
electrons was above the noise level. The contamination 
of the pion sample by other hadrons as well as leptons is 
negligible. 

The RICH detector identifies pions, kaons, and protons 
in the momentum range 2 GeV < p < 15 GeV. In the 
semi-inclusive analysis reported in this paper a momen- 
tum range of 4 GeV < p < 13.8 GeV was used for consis- 



exp 



qthco;a,5 



(27) 



a particle hypothesis i — t:,K,p with the largest total 
likehhood £*°* = Cf-£f is assigned to each hadron track. 

Identification efhciencies and probabilities for contam- 
ination of hadron populations from misidentification of 
other hadrons were estimated with a Monte Carlo sim- 
ulation which had been calibrated with pion and kaon 
tracks from experimentally reconstructed (j), me- 
son and A hyperon decays. In the analysis, each pion and 
kaon track was assigned a weight cof'^ accordingly. The 



15 



number of counts of pions and kaons, 

i 

were computed as the sums of these weights. 

V. ASYMMETRIES 
A. Measured Asymmetries 



In the data sample that satisfied the data quahty cri- 
teria described in the previous section, events were se- 
lected for analysis if they passed the DIS trigger (see sec- 
tion ^Oj. Tracks with a minimum energy of 3.5 GeV in 
the calorimeter that were identified as leptons by the PID 
system were selected as candidates for the scattered DIS 
particle by imposing additional requirements on the track 
kinematics. A requirement of > 1 GeV^ selected hard 
scattering events. Events from the nucleon resonance re- 
gion were eliminated by requiring > lOGeV^. The 
requirement y < 0.85 reduced the number of events with 
large radiative corrections. 

Hadron tracks coincident with the DIS positron were 
identified as semi-inclusive hadrons if the fractional en- 
ergy z of the hadron was larger than 0.2 and xp was 
larger than 0.1. These limits suppress contributions from 
target fragmentation. Hadrons from exclusive processes, 
such as diffractive vector meson production, were sup- 
pressed by requiring that z be smaller than 0.8. 

The final statistics of inclusive DIS events and SIDIS 
hadrons are given in Tab. II VI The numbers are presented 
in terms of statistically equivalent numbers of events N^q. 
This quantity is the number of unweighted events with 
the same relative error as the sum of weighted events iV, 



eq 



eq 



(29) 



The weights uji are defined in section IIVDI for hadrons 
identified by the RICH detector in semi-inclusive events. 
For pions identified by the threshold Cerenkov counter, 
for undifferentiated hadrons, and for inclusive DIS events 
oji = 1. An additional weight factor of ±1 is applied 
according to the event classification as signal or charge- 
symmetric background (see below). 



TABLE IV: Statistically equivalent number of counts of DIS 
events and SIDIS hadrons for the hydrogen and the deuterium 
data. 









SIDIS events 




Target 


DIS evts. 


7V+ 


TT" K+ 


K- 


H 


1.7 X 10^ 


117 X 10^ 


82 X 10^ 




D 


6.7 X 10® 


491 X 10^ 


385 X 10^ 76 X 10^ 


33 X 10^ 



The inclusive (semi-inclusive) data samples were used 
to calculate the measured positron-nucleon asymmetry 



(28) ^ii''''' in bins of x (or z). 



A 



(h) 



(h) 



(h) 



(30) 



Here {^^) is the number of DIS events for tar- 
get spin orientation parallel (anti-parallel) to the beam 

spin orientation, and Nf^ i-^iT) ^he corresponding 
numbers of semi-inclusive DIS hadrons. The luminos- 
ity (L^) for the parallel (anti-parallel) spin state is 

corrected for dead time, while (Lp) is the luminos- 
ity corrected for dead time and weighted by the product 
of beam and target polarizations for the parallel (anti- 
parallel) spin state. Values for the beam and target po- 
larizations are given in Tabs. ^ and IIIII The bins in x 
used in the analysis are defined in Tab. 



TABLE V: The bins in x used in the analyses presented in 
this paper. 



bin 


1 


2 


3 


4 


5 


6 


7 


8 


9 


^low 
3^up 


0.023 
0.040 


0.040 
0.055 


0.055 
0.075 


0.075 
0.1 


0.1 
0.14 


0.14 
0.2 


0.2 
0.3 


0.3 
0.4 


0.4 
0.6 



In the deuteron — a spin-1 particle — another polar- 
ized structure function bf arises from binding effects as- 
sociated with the D-wave component of the ground state 
|43| . This structure function may contribute to the cross 
section if the target is polarized with a population of 
states with spin-projection Sz — that is not precisely 
1/3, i.e. a substantial tensor polarization. Because the 
maximum vector polarization can only be accomplished 
with a high tensor polarization in a spin 1 target, mea- 
surements in HERMES, of necessity, can include signif- 
icant contributions from the tensor analyzing power of 
the target. For inclusive scattering, the spin asymmetry 
is of the form 

1, 



0'l/2 - <^3/2 I 



-TAt] 



(31) 



0'l/2 + ^3/2 

where T is the tensor polarization, and At is the tensor 
analyzing power. The bi structure function is measured 
by At, i.e. At « 26i/3Fi. Studies by the HERMES 
collaboration indicate that bf is small |4J|, and that the 
tensor contribution to the the inclusive deuteron asym- 
metry is less than « 0.5 — 1.0% of the measured asymme- 
try. For this reason, tensor contributions were assumed 
to be negligible for all the spin asymmetries presented 
here. 



B. Charge— Symmetric Background 



The particle count rates were corrected for charge sym- 
metric background processes {e.g. j —^ e'^ + e~). The 
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rate for this background was estimated by considering 
lepton tracks with a charge opposite to the beam charge 
that passed the DIS restrictions. It was assumed that 
these leptons stemmed from pair-production processes. 
The rate for the charge symmetric background process 
(where the particle is detected with the same charge as 
the beam but originating from pair production) is the 
same. The number of events with an opposite sign lep- 
ton is therefore an estimate of the number of charge sym- 
metric events that masquerade as DIS events. They were 
subtracted from the inclusive DIS count rate. Hadrons 
that were coincident with the background DIS track and 
that passed the SIDIS limits were also subtracted from 
the corresponding SIDIS hadron sample. The DIS back- 
ground rate was ~ 6 % with respect to the total DIS rate 
in the smallest x-bin, falling off quickly with increasing 
X. The overall background fraction from this source was 
1.4%. 



C. Azimuthal Acceptance Correction 

The measured semi-inclusive asymmetries were cor- 
rected for acceptance effects due to the nonisotropic az- 
imuthal acceptance of the spectrometer. These accep- 
tance effects arise because of an azimuthal dependence of 
the polarized and unpolarized semi-inclusive cross sec- 
tions due, e.g . to non-zero intrinsic transverse parton 
momenta |4!tI |. Taking into account the azimuthal de- 
pendence, the measured semi-inclusive asymmetry given 
in Eq. (jSOJl is modified [i^ . 



A'f-f gi(cos( 



ILL 



1 + g^(cos( 



/uu 



(32) 



where {cos 4>)ll and {coa 4))uu are the cos0 moments of 
the semi-inclusive polarized and unpolarized cross sec- 
tion, respectively and Co and Ci are the lowest order 
Fourier coefficients of the spectrometer's azimuthal ac- 
ceptance. 

The semi-inclusive asymmetries were corrected for 
the unpolarized moment (cos (j))uu using a determination 
of the moments from HERMES data. The polarized mo- 
ment (cos </))ll was found to be neghgible as expected in 
Ref. j4g. The correction for the unpolarized moment to 
the asymmetries is 10% for x < 0.1. In the measured x 
range the absolute correction of the semi-inclusive asym- 
metries is small, because of the small size of the asym- 
metries at low X (see below) and because the correction 
is small for x > 0.1. The correction to the asymmetries 
as function of z is about 10% at small z and becomes 
smaller for larger values of z. 



D. Radiative and Detector Smearing Effects 

The asymmetries were corrected for detector smearing 
and QED radiative effects to obtain the Born asymme- 



tries which correspond to pure single photon exchange 
in the scattering process. The corrections were applied 
using an unfolding algorithm that accounts for the kine- 
matic migration of the events. As opposed to iterative 
techniques described in Ref. {47| for example, this al- 
gorithm does not require a fit of the data. The final 
Born asymmetries shown in Figs . [T5I and 1 1 41 depend only 
on the measured data, on the detector model, on the 
known unpolarized cross sections, and on the models for 
the background processes. Another advantage is the un- 
ambiguous determination of the statistical variances and 
covariances on the Born asymmetries based on the sim- 
ulated event migration. A description of the unfolding 
algorithm and the input Monte Carlo data is found in 
App.El 




FIG. 11: The observed asymmetries Ay and AJ^ on the pro- 
ton target compared with the corresponding Born asymme- 
tries. The Born asymmetries are offset horizontally for better 
presentation. See the text for details. 
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FIG. 12: Uncertainty inflation caused by detector smearing 
and QED radiative effects. The left hand panel shows the un- 
certainty inflation of the inclusive asymmetry for the proton 
and the right hand panel that of the positive pion asymme- 
try. In both panels the open triangles present the uncertainty 
inflation caused by detector effects, the open squares present 
the inflation caused by QED radiative effects, and the fllled 
circles show the total uncertainty inflation. 
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FIG. 13: The inclusive and semi-inclusive Born level asymmetries on the proton, corrected for instrumental smearing and QED 
radiative effects. The error bars give the statistical uncertainties, and the shaded bands indicate the systematic uncertainty. 
The open squares show the positive and negative hadron asymmetries measured by the SMC collaboration, limited to the 
HERMES i-range 



The impact of the unfolding on the asymmetries is il- 
lustrated with the inclusive and the positive pion asym- 
metries on the proton in Fig. 1111 The unfolding pro- 
cedure shifts the central asymmetry values only by a 
small amount. This is expected for cross section asym- 
metries. Smearing results in a loss of information about 
more rapid fluctuations that may be present in the data. 
Therefore correcting for this loss by effectively enhancing 
"higher frequency" components inevitably results in an 
inflation of the uncertainty of each data point. The un- 
certainty inflation introduced by the unfolding is shown 
in Fig. 1121 The uncertainty at low x is significantly in- 
creased by the QED background in the case of the in- 
clusive asymmetry. At large values of x and in the case 
of the semi-inclusive asymmetries, the inflation by QED 
radiation is mainly due to inter-bin event migration. De- 
tector smearing effects, which are largest at high x, in- 
crease the uncertainties through inter-bin migration and 
a small number of events that migrate into the accep- 
tance. Uncertainty inflation due to interbin migration 
increases rapidly as the bin size is reduced to be compa- 
rable to the instrumental resolution. 

It should be noted that this unfolding procedure is 
more rigorous than the procedure applied in previous 
DIS experiments and previous analyses of this experi- 
ment. The size of the uncertainties is larger in the cur- 
rent analysis due to the explicit inclusion of correlations 
between x-bins and the model-independence of the un- 



folding procedure. This should be borne in mind when 
comparing the current data to other results. 



E. Results for the Asymmetries 

The asymmetries are related to the inclusive 

and semi-inclusive photon-nucleon asymmetries A^''' 
through the kinematical factors 77 and 7 and the depolar- 
ization factor D (see Eq. (|13ll '). The Born level asymme- 
tries A^'''' on the proton and on the deuteron targets are 
shown in Figs. 1131 and 1141 and are listed in Tabs. IXITI and 
IXIIII in Add. IbI respectively. The present results on the 
proton target supersede earher results pubhshed in [i^ . 

The inclusive asymmetries on the proton and the 
deuteron are determined with high precision. On both 
targets, they are large and positive. A detailed discus- 
sion and a determination of the spin structure functions 
gi from inclusive scattering data is given in Ref. 50] for 
the proton and a forthcoming paper in the case of the 
deuteron. 

The asymmetries for undifferentiated positive and neg- 
ative hadrons on both targets are compared with mea- 
surements performed by the SMC collaboration . The 
statistical uncertainties of the HERMES data are signif- 
icantly better than those of the SMC data. Pion asym- 
metries on the proton and pion and kaon asymmetries 
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FIG. 14: The inclusive and semi-inclusive Born level asymmetries on the deuteron. One data point at a; = 0.45 for the K~ 
asymmetry including its large error bar is outside the displayed range; all data points are listed in Tab. IXIlTI See Fig. Il.'-il for 
details. 



on the deuteron were measured for the first time. The 
pion asymmetries are determined with good precision, 
whereas the kaon asymmetries have larger statistical un- 
certainties. Except for the K~ asymmetry, all asymme- 
tries are seen to be mostly positive, which is attributed 
to the dominance of scattering off the u-quark. The frag- 
mentation into negative kaons (us-mesons) has in com- 
parison to the other hadrons an increased sensitivity to 
scattering off u and s-quarks, which makes the K~ asym- 
metry a useful tool to determine the polarization of these 
flavors. 



F. z— Dependence of the Asymmetries 

Because the ratio of favored to unfavored fragmenta- 
tion functions is known to vary substantially with z, a 
z-dependence of the asymmetries could be induced by 
the variation of the relative contributions of the various 
quark flavors to fragmentation. The observation of a z- 
dependence of the asymmetries could also be caused by 
hadrons in the semi-inclusive data sample that originate 
from target fragmentation as opposed to current frag- 
mentation, which is associated with the struck quark. 
Furthermore, hadrons from non-partonic processes such 
as diffractive interactions could play an important role in 
the semi-inclusive DIS data sample 51]. For example, at 
high fractional energies z, it is possible that hadrons from 
exclusive processes are misinterpreted as SIDIS hadrons. 
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FIG. 15: The semi-inclusive Born asymmetries for positive 
and negative pion production on the proton as a function of 
z. The error bars indicate the statistical uncertainties and the 
error band represents the systematic uncertainties. The solid 
line is the z dependence from the Monte Carlo simulation of 
the asymmetries. 



To explore these possibilities, and to test the Jet- 
SET fragmentation model used here (see section IVI A|) 
in the Monte Carlo simulation of the scattering process, 
the semi-inclusive asymmetries were extracted in bins 
of z. They were calculated with the same kinematical 
limits described above, except for the requirement on 
xp, which is highly correlated with the limit on z and 
was therefore discarded. Events were accepted over the 
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range 0.023 < x < 0.6. The semi-inclusive pion asym- 
metries for the proton are shown in Fig. 1151 together with 
a curve of the asymmetries from the Monte Carlo simu- 
lation. The agreement between experimental and simu- 
lated data provide confirmation that the fragmentation 
process is consistently modeled. 



G. Systematic Uncertainties in Ai 

Systematic uncertainties in the observed lepton- 
nucleon asymmetries a'^^'^ arise from the systematic un- 
certainties in the beam and target polarizations. The un- 
folding of the observed asymmetries also increases these 
uncertainties. A systematic uncertainty due to the RICH 
hadron identification was estimated to be small as the ef- 
fect of neglecting the hadron misidentification (neglecting 
the off-diagonal elements of w appearing in Eq. H28|) ) was 
found to be negligible. Therefore, it was not included in 
the semi-inclusive deuterium asymmetries. 



TABLE VI: The fractional systematic uncertainties on A\ 
averaged over x. 



Source 


Hydrogen data 


Deuterium data 


Beam polarization 


4.2% 


2.3% 


Target polarization 


5.1% 


5.2% 


Azimuthal acc. (SIDIS) 


3.0% 


3.1% 


QED rad. corr. (DIS) 


2.0% 


2.0% 


QED rad. corr. (SIDIS) 


1.0% 


1.0% 


Detector smearing 


2.0% 


2.0% 


R 


1.1% 


1.1% 


92 


0.6% 


1.4% 



Additional uncertainties arise due to the finite MC 
statistics, when the corrections for detector smearing and 
QED radiation are applied. They are included in the 
statistical error bars in the figures and are listed in a 
separate column in the tables shown in the appendix. 

In forming the photon-nucleon asymmetries A^"^ , sys- 
tematic uncertainties due to the parameterization of the 
ratio R and the neglect of the contribution from the 
second polarized structure function g2 were included 
[43. IH^ . The relative systematic uncertainties are sum- 
marized in Tab. IVII The total systematic uncertainties 
on the asymmetries are shown as the error bands in the 
figures. 

The interpretation of the extracted asymmetries may 
be complicated by contributions of pseudo-scalar mesons 
from the decay of exclusively produced vector mesons, 
mostly p^'s producing charged pions. The geometric ac- 
ceptance of the spectrometer is insufficient to identify and 
separate these events, as typically only one of the decay 
mesons is detected. However, the fractional contribu- 
tions of diffractive vector mesons to the semi-inclusive 
yields were estimated using a Pythia6 event genera- 



tor JSJ that has been tuned for the HERMES kinemat- 
ics |5J|. The results range from 2%(3%) at large x to 
10%(6%) at small x for pions(kaons) for both proton and 
dcuteron targets. Although some data of limited preci- 
sion for double-spin asymmetries in and (j> production 
have been measured by HERMES jS^ , no information is 
available on the effects of target polarization on the angu- 
lar distributions for the production and decay of vector 
mesons. Therefore it was not possible at this time to 
estimate the effect of the decay of exclusively produced 
vector mesons on the semi-inclusive asymmetries. 

The measurement of asymmetries as opposed to total 
cross sections has the advantage that acceptance effects 
largely cancel. Nevertheless, the forward acceptance of 
the spectrometer restricts the topology of the DIS elec- 
tron and the SIDIS hadron in the final state. It was sug- 
gested ISflj that a resulting cutoff in transverse hadron 
momentum leads to a bias in the contributions of photon 
gluon fusion (PGF) and QCD Compton (QCDC) pro- 
cesses to the total DIS cross section. This bias could lead 
to an incorrect measurement of the polarizations of the 
quarks using SIDIS asymmetries. The momentum cut 
(4GeV < p < 13.8 GeV) on the coincident hadron tracks 
for particle identification using the Cerenkov/ RICH (cf. 
section llV Dp could potentially introduce further bias. 
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FIG. 16: Born level Monte Carlo asymmetries on the proton 
in the experimental acceptance and in An. The left hand 
plot compares the semi-inclusive asymmetries AJ^p and the 

right hand plot the semi-inclusive asymmetries AJ^p . The 
asymmetries in the experimental acceptance also include the 
hadron momentum cut. For display purposes the full data 
points have been offset horizontally. The lower panels present 
the same data in the form of the difference in the asymmetries 
divided by the total experimental uncertainty ai,orn in the 
corresponding measured Born asymmetry. 



Possible effects on the asymmetries due to the accep- 
tance of the HERMES spectrometer were studied with 
the HERMES Monte Carlo simulation. Born level data 
were generated using a scenario in which contributions 
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to the cross sections from PGF and QCDC processes 
were found to be smaller than 7% and 18% respec- 
tively. These values were obtained in a scheme of cut 
offs against divergences in the corresponding QCD ma- 
trix elements which require quark-antiquark pairs to have 
masses iriq-q > 1 GeV and > 0.005 W^. This is to be com- 
pared with the default values used in the purity analy- 
sis of niqq > 2 GeV and > 0.005 W^^. In this default 
case the contributions from PGF and QCDC processes 
were less than 1.5% and 3% respectively. In the sce- 
nario employed for the acceptance study the effect of the 
experimental acceptance was determined to be negligible 
compared to the uncertainties in the data. The semi- 
inclusive 7r+ and tt^ asymmetries in An and inside the 
acceptance are compared in Fig. ^1 Acceptance effects 
on the Born level asymmetries are small and corrections 
are not necessary. 



VI. QUARK HELICITY DISTRIBUTIONS 

A. Quark Polarizations and Quark Helicity 
Densities 

A "leading order" analysis which included the PDF 
and QCDC processes discussed in the previous section 
was used to compute quark polarizations from the Born 
asymmetries. The contribution of exclusively produced 
vector mesons is not distinguished in this extraction. The 
analysis based on Eq. (|19|l combines the Born asymme- 
tries in an over-constrained system of equations, 



on the deuteron: 



(33) 



where the elements of the vector Ai{x) are the measured 
inclusive and semi-inclusive Born asymmetries and the 
vector Q{x) contains the unknown quark polarizations. 
The matrix TV is the nuclear mixing matrix that accounts 
for the probabilities for scattering off a given nucleon in 
the deuteron nucleus and the nucleon's relative polariza- 
tion. The matrix V contains as elements the effective 
purities for the proton and the neutron. These elements 
were obtained by integrating Eq. (|19|l over the range in 
z and giving 



Ag(a;) 
q{x) 



(34) 



where Vg{x) is now the effective spin-independent purity 



(35) 



The vector Ai (x) includes the inclusive and the semi- 
inclusive pion asymmetries on the proton, and the inclu- 
sive and the semi-inclusive pion and kaon asymmetries 



Ai{x) = (Ai,p(a;), Afpix), 



(36) 



The semi-inclusive asymmetries of undifferentiated 
hadrons were not included in the fit because they are 
largely redundant with the pion and kaon asymmetries 
and thus do not improve the precision of the results. 

The nuclear mixing matrix Af combines the proton and 
neutron purities into effective proton and deuteron puri- 
ties. The relation is trivial for the proton. In the case of 
the deuteron purities, the matrix takes into account the 
different probabilities for scattering off the proton and 
the neutron as well as the effective polarizations of the 
nucleons. The probabilities were computed with hadron 
multiplicities measured at HERMES and using the NMC 
parameterization of F2 .57.] . The D-state admixture in 
the deuteron wave function of (5 ± 1) % js^,!!^ leads to 
effective polarizations of the nucleons in the deuteron of 
Pp,D = Pn,D = (0.925 ± 0.015) PD. 

The purities depend on the unpolarized quark densi- 
ties and the fragmentation functions. The former have 
been measured with high precision in a large number 
of unpolarized DIS experiments. The CTEQ5L parton 
distributions ,60,J incorporating these data were used in 
the purity determination. Much is known about frag- 
mentation to mesons at collider energies. However, the 
application of this information to fixed target energies 
presents difficulties, especially regarding strange frag- 
mentation, which at lower energies no longer resembles 
that of lighter quarks. A recent treatment 61] using ex- 
tracted fragmentation functions in an analysis of inclu- 
sive and semi-inclusive pion asymmetries for the proton 
demonstrates the shortcomings of using the limited ava- 
iable data base for fragmentation functions. Hence the 
interpretaton of the present asymmetry data requires a 
description of fragmentation that is constrained by me- 
son multiplicities measured at a similar energy. Such 
multiplicities within the HERMES acceptance are avail- 
able, but those for kaons not yet available corrected to 
47r acceptance. Hence the approach taken here was to 
tune the parameters of the LUND string model imple- 
mented in the Jetset 7.4 package (62] to fit HERMES 
multiplicities as observed in the detector acceptance. 

In the LUND model, mesons are generated as the 
string connecting the diquark remnant and the struck 
quark is stretched. Quark-antiquark pairs are generated 
at each breaking of the string. Even though the leading 
hadron is often generated at one of the string breaks and 
not at the end, the flavor composition of any hadron ob- 
served at substantial z retains a strong correlation with 
the flavor of the struck quark. It is this correlation which 
provides semi-inclusive flavor tagging. Contrary to some 
speculations based on a misunderstanding of the Jetset 
code [i^ , this feature of the Lund string model is inde- 
pendent of in lepton-nucleon scattering. The quarks 
associated with either half of the string retain the infor- 
mation on the flavor of the struck quark. The LUND 
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FIG. 17: Multiplicities of charged pions and kaons in the 
HERMES acceptance compared with Monte Carlo data as 
function of z. Statistical uncertainties in the data are too 
small to be visible. The multiplicities shown were not cor- 
rected for radiative and instrumental effects. 



model has proved to be a reliable widely accepted means 
of describing the fragmentation process. 

The string breaking parameters of the LUND model 
were tuned to fit the hadron multiplicities measured at 
HERMES in order to achieve a description of the frag- 
mentation process at HERMES energies . A compar- 
ison of the measured and the simulated hadron multiplic- 
ities is shown in Fig. El The tuned Monte Carlo sim- 
ulation reproduces the positive and negative pion mul- 
tiplicities and the negative kaon multiplicities while the 
simulated positive kaon multiplicities are smaller than 
those measured. A similar disagreement is also reported 
by the EMC experiment [Hi. 

The purities were computed in each x-bin i from the 
described tuned Monte Carlo simulation of the entire 
scattering process as 



(37) 



In this expression, iV^ is the number of hadrons of type h 
in bin i passing all kinematic restrictions when a quark of 
flavor q was struck in the scattering process. The purities 
include effects from the acceptance of the spectrometer. 
In Fig. llSl the purities for the proton and the neutron are 
shown. It is evident from these plots that the w-quark 
dominates the production of hadrons, due to its charge 
of 2/3 and its large number density u{x) in the nucleon. 



In particular the large contribution by the it-quark to 
7r+ production from both proton and neutron targets 
provides excellent sensitivity to the polarization of the 
M-quark. The d-quark becomes accessible through the 
production of negative pions, which also separates the u 
and d flavors. More generally, contributions from the sea 
quarks can be separated from the valence quarks through 
the charge of the final state hadrons. Finally the mea- 
surement of negative kaons is sensitive to strange quarks 
and the anti-strange quark can be accessed through posi- 
tive kaons in the final state. However, large uncertainties 
in the strange sea distributions are expected, because the 
strange and anti-strange purities are small in comparison 
to those for the other flavors. Some of the strange quark 
purities vary rapidly with Q^. 

The quark polarizations [Aq/q](x) are obtained by 
solving Eq. H33() . A combined fit was carried out for all 
x-bins to account for the statistical correlations of the 
Born asymmetries (cf. Eq. ljA7p ). Accordingly, the vec- 
tors Ai and MVQ are arranged such that they contain 
consecutively their respective values in all a;-bins, 



{Ai{xi), Ai{x2), Ai{xq)^ 



(38) 



and analogously for NVQ. The polarizations follow by 
minimizing 



X 



= [Ai-NVQ] Va'[Ai-NVQ 



(39) 



where Va is the statistical covariance matrix (Eq. ljA7|l ') 
of the asymmetry vector Ai. It accounts for the corre- 
lations of the various inclusive and semi-inclusive asym- 
metries as well as the inter-bin correlations. 

The systematic uncertainties of the asymmetries were 
not included in the calculation of . The dominant con- 
tribution to these uncertainties arises from the beam and 
target polarizations, which affect the asymmetries in a 
nonlinear manner. It is natural to linearly approximate 
these contributions as off-diagonal inter-bin correlations 
in the systematic covariance matrix of the set of asym- 
metries for all bins. However, when such a matrix is 
included in the fit based on linear recursion, the inac- 
curacies in this linearization were found to introduce a 
significant bias in the fit. Hence, the systematic uncer- 
tainties were excluded from the fit, but were included in 
the propagation of all of the uncertainties in the asym- 
metries into those on the results of the fit. 

It was found that the data do not significantly con- 
strain As(a;). The results presented here were extracted 
with the constraint As(a;) = 0. A comparison of the 
fit using this constraint with a fit without assumptions 
on the polarizations of the quark flavors showed that 
the constraint had negligible impact on the final re- 
sults for the unconstrained flavors and their uncertain- 
ties. In addition the resulting polarizations were found 
to be in good agreement with the results of a fit under 
the assumption of a symmetrically polarized strange sea 
[As/s](x) = [As/s](x). 




FIG. 18: Purities for positive and negative pions and kaons on a proton and a free neutron target. Each column corresponds to 
scattering off a certain quark flavor. The shaded bands indicate the estimated systematic uncertainties due to the fragmentation 
parameters. The estimate derives from a comparison of two parameter sets. See text for details. 



Assuming an unpolarized anti-strange sea the vector 
of polarizations Q{x) in each x-bin is given by 



Qix) 



Au. . Ad. . Au, . Ad. . As, , 

(x), -r[x), —[x), -r(a;), — (x) 

u a u d s 



(40) 

As a further constraint the polarizations of the u, d, and 
s-quarks were fixed at zero for values of x > 0.3. The 
effects of this and of fixing the s polarization at zero were 
included in the systematic error. The constraints reduced 
the number of free parameters by fifteen, leaving 39 pa- 
rameters in the fit. The solution obtained by applying 
linear regression is 

Q = irIf{VA)-'Vefy'vJf{VA)-'A[, (41) 

where A[ = Ai -NVQf^x, Vef = J^V, and Qfi^ is 
the set of constrained polarizations. The covariance ma- 
trix of the quark polarizations propagated from the Born 
asymmetries is 



v(Q)= {piA^Ar'v.fy'pfyvAY 



(42) 



where the covariance matrix V^°* includes the statisti- 
cal and the systematic covariances, Vjf* = Va + . 
The resulting solution is shown in Fig. El The value of 
the x^/ndf of the fit is 0.91. The reasonable value 
confirms the consistency of the data set with the quark 
parton model formalism of section II. C. Removing the 
inclusive asymmetries from the fit has only a small effect 
on the quark polarizations and their uncertainties. 

The polarization of the u-quarks is positive in the mea- 
sured range of x with the largest polarizations at high x 
where the valence quarks dominate. The polarization of 
the d-quark is negative and also reaches the largest (neg- 
ative) polarizations in the range where the valence quarks 
dominate. The polarization of the light sea flavors u and 
d, and the polarization of the strange sea are consistent 
with zero. The values of jndj for the zero hypotheses 
are 7.4/7, 11.2/7, and 4.3/7 for the m, the d, and the 
s-quark, respectively. 

The quark polarizations in Fig. E| are presented at 
the measured Q^-values in each bin of x. The Q^- 
dependence is predicted by QCD to be weak and the 
inclusive and semi-inclusive asymmetries measured by 
HERMES (cf. Figs, ini m and Ref. ^49,]) and SMC Hi 
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FIG. 19: The quark polarizations in the 5 parameter x 9 x- 
bins fit. The polarizations, shown as a function of x, were 
computed from the HERMES inclusive and semi-inclusive 
asymmetries. The error bars are the statistical uncertain- 
ties. The band represents the total systematic uncertainty, 
where the light gray area is the systematic error due to the 
uncertainties on the fragmentation model, and the darlc gray 
area is from the contribution of the Born asymmetries. 



at very different average show no significant Q^- 
dependence when compared to each other. The quark 
polarizations [Aq/q]{x) are thus assumed to be Q^- 
independent. 

The quark hehcity densities Aq{x, Qq) are evaluated 
at a common Qq = 2.5 GeV^ using the CTEQ5L un- 
polarized parton distributions. Because the CTEQ5L 
compilation is based on fits to experimental data for 
^2(2;), the relationship between F2{x) and Fi{x) as given 
by Eq. is here taken into account. The factor 
Cr = (1 + i?)/(l + 7^) connects CTEQ5L tabulations 
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FIG. 20: The quark helicity distributions x Aq{x,QQ) eval- 
uated at a common value of Qg = 2.5 GeV^ as a function of 
X. The dashed line is the GRSV2000 parameterization (LO, 
valence scenario) scaled with 1/(1 + R) and the dashed- 
dotted hue is the Bliimlein-Bottcher (BB) parameterization 
(LO, scenario 1) [6^ . See Fig. 1191 for explanations of the 
uncertainties shown. 



with the parton distributions q{x) required here. In the 
present analysis the parameterization for R{x, Q^) given 
in Ref. [G^I was used. The results are presented in Fig.EUl 
The data are compared with two parton helicity distri- 
butions 18, "66] derived from LO fits to inclusive data. 
The GRSV2000 parameterization, which was fitted using 
the assumption i? = 0, is shown with the scaling factor 
1/(1 + i?) to match the present analysis. While in the 
Bliimlein-Bottcher (BB) analysis equal helicity densities 
for all sea flavors are assumed, in the GRSV2000 "va- 
lence fit" a different assumption is used, which leads to 
a breaking of fiavor symmetry for the sea quark helicity 
densities. In Tab. lVIll the x^-values of the comparison of 
the measured densities with these parameterizations and 
the zero hypothesis are given. The measured densities 
are in good agreement with the parameterizations. The 
data slightly favor the BB parameterization of the u~ and 
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TABLE VII: Comparison of the measured quark helicity den- 
sities with the parameterizations obtained from LO QCD fits 
to inclusive data and with the zero hypothesis. Listed are the 
reduced values /ndf for each hypothesis. 





xAu{x) 


xAd{x) 


x/S.u{x) 


a;Aj(x) 


xl^s{x] 


GRSV2000 val. 


1.45 


0.93 


1.54 


1.44 


0.60 


BB (scenario 1) 


1.02 


1.06 


0.97 


1.32 


0.95 


x/S.q{x) = 


13.19 


2.50 


1.06 


1.60 


0.61 


ndf 


9 


9 


7 


7 


7 



M-flavors, while for the other flavors the agreement with 
both parameterizations is equally good. Within its un- 
certainties the measured strange density is in agreement 
with the very small non-zero values of the parameteriza- 
tions as well as with the zero hypothesis. 

The total systematic uncertainties in the quark po- 
larizations and the quark helicity densities include con- 
tributions from the input asymmetries and systematic 
uncertainties on the purities, which may arise from the 
unpolarized parton distributions and the fragmentation 
model. Since the apphed CTEQ5L PDFs [63| are pro- 
vided without uncertainties, no systematic uncertainty 
from this source was assigned to the purities. 

The uncertainties of the fragmentation model would 
be ideally calculated by surveying the (unknown) x^" 
surface of the space of Jetset parameters that were used 
to tune the Monte Carlo simulation. At the time of pub- 
lication such a computationally intensive scan was not 
available. Instead the uncertainties were estimated by 
comparing the purities obtained using the best tune of 
Jetset parameters described above to a parameter set 
which was derived earlier [49l | . This earlier parameter set 
was also obtained from a similar procedure of optimizing 
the agreement between simulated and measured hadron 
multiplicities. However, because of the lack of hadron 
discrimation in a wide momentum range before the avail- 
ability of the RICH detector and limited available com- 
puter power, this earlier parameter tune optimized only 
three Jetset parameters, while in the current tune 
eleven parameters were optimized from their default set- 
tings. The differences in the resulting purities from using 
these two different tunes of Jetset parameters are shown 
as the shaded bands in Fig. 1181 

The contributions from this systematic uncertainty es- 
timate on the purities to the total systematic uncertain- 
ties of the resulting helicity densities and quark polar- 
izations are shown as the light shaded bands in Figs. 
1201 and Fig. |^ In the case of the u and d quark, the 
resulting uncertainty contributions due to the fragmen- 
tation model are small compared with those related to 
the systematic uncertainties on the asymmetries. They 
are of equal or larger size in the case of the sea quarks 
and they dominate in case of the light sea quark helicity 
difference discussed in the following subsection. 



B. Isospin Asymmetry in Helicity Densities of 
Light Sea Quarks 

In the unpolarized sector, the breaking of flavor sym- 
metry for the light sea quarks {d{x) — u{x) > 0) and 
consequently the violation of the Gottfried sum rule is ex- 
perimentally well established [HjH^It^ and is described 
by various non-perturbative models. Two comprehensive 
reviews of these models can be found in 71, 72] . Such 
models also predict a flavor asymmetry in the light sea 
helicity densities. Sizeable asymmetries of similar mag- 
nitude arepredicted by the chiral quark soliton model 
(xQSM) IH El El Izi El, m, which is an effective 
theory where baryons appear as soliton solutions of the 
chiral Lagrangian, a statistical model E3 ^^at describes 
the nucleon as a gas of masslcss partons, models based 
on the meson cloud picture ,8^ ,81,] as well as the chiral 
chromodielectric model, which is a bag-like confinement 
model 's^. The meson cloud model published in 81] de- 
viates most from all other mentioned models as it predicts 
an asymmetry which is smaller in magnitude but has the 
opposite sign Au{x) — Ad{x) < to the other models, 
which all predict a positive value of this quantity. 

The measured (semi-)inclusive asymmetries discussed 
above were used in a modified fit to compute this flavor 
asymmetry. In this fit the parameter [AM/?2](a;) was re- 
placed by [(Am — A(J) / (u — d)] (x) and the system of hnear 
equations H33|l was solved for the following vector Q'{x) 
of quark polarizations: 

<3'(i) = 

\ u d u — d d s J 

(43) 

The fiavor asymmetry x [Au(x) ~ Ad(2:)] is presented in 
Fig. 1^ For c omp arison the same quantities calculated 
m the xQSM III and in a meson cloud model ISlI are 
shown. These models were chosen for presentation as 
they indicate the most positive and most negative pre- 
dictions. All other cited model predictions yield similar 
x-dependences and positive values Au(a;) — Ad[x) > 0. 
Because of the close similarity of these model curves when 
plotted at the scale of Fig. [2] they are not displayed in 
this figure. 

The value of the 'x^ /ndf for the symmetry hypothe- 
sis A?2(a;) = Ad{x) is 7.7/7. The /ndf values of the 
comparisons with the predictions shown for the xQSM 
and for the meson cloud model are 17.6/7 and 8.1/7, 
respectively, where the value for the xQSM takes into 
account its uncertainty. This analysis of the HERMES 
data therefore favors a symmetric polarized light flavor 
sea and exclude the prediction of the xQSM at the 97% 
confldence level. 
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FIG. 21: The flavor asymmetry in the helicity densities of 
the light sea evaluated at Qq — 2.5 GeV^. The data are com- 
pared with predictions in the xQSM [t^ and a meson cloud 
model 81]. The solid line with the surrounding shaded band 
show the xQSM prediction together with its ±1 a uncertain- 
ties while the dash-dotted line shows the prediction in the 
meson cloud model. The uncertainties in the data are pre- 
sented as in Fig. 1191 



VII. As FROM THE ISOSCALAR 
EXTRACTION METHOD 

The strange quark polarization in the nucleon is one 
of the most interesting quantities that can be determined 
from SIDIS data. Kaon asymmetries provide the largest 
sensitivity to the strange quark polarization because the 
kaon contains a valence strange (anti-)quark. Unfortu- 
nately, experimental data on separate favored and dis- 
favored fragmentation of strange quarks into kaons are 
scarce at best in the kinematic region of HERMES, 
thereby introducing large systematic uncertainties on the 
strange quark purities (cf. Fig. [TH|l . This section presents 
an alternative approach for the extraction of the strange 
quark polarization, which uses only data on the asym- 
metry of the total {K^ + K^) charged kaon flux for a 
deuterium target, and which does not rely on a Monte 
Carlo model of the fragmentation process. 

The total strange quark helicity density AS'(x) = 
As(x) -|- As(x) carries no isospin. It can hence be ex- 
tracted from scattering data off a deuterium target alone, 
which is isoscalar. Furthermore, the fragmentation func- 
tions -D^jj^J^ for the total kaon flux were measured at 
e^e~ collider experiments with satisfactory precision. 

The analysis is performed as a two component ex- 
traction of the non-strange and total strange quark po- 
larizations, [l^Q/Q]{x) and [AS'/5](a;). Here, Q{x) = 
u{x) -\- u{x) -f d{x) -f d{x). Only two measured asymme- 
tries, the inclusive Ai d and the semi-inclusive 
asymmetries on the deuterium target, are used for the 
extraction. The same purity approach detailed in sec- 



tions |n| and IVII is used, except in the present case the 
purity matrix 'P{x) consists of a 2 x 2 matrix only: 



V{x) = 



P, 



Pq{x) Psix) 

K+ + K- /\ r,K+ + K- 



Q 



{x) 



(44) 



Due to row-wise unitarity, there are only two indepen- 
dent elements in this matrix. 

Based on the wealth of e+e" collider data, parameteri- 
zations of the total charged kaon fragmentation functions 

P^^ jt^ '5^) are available as a function of z and 
|8a. |8J| . The light quarks u, d, s are assumed to be mass- 
less in the analyses of the e+e" colHder data |8^ 8^. 
This is justifled in view of the high center-of-mass en- 
ergies available in these experiments. However, this ap- 
proximation is not valid for the lower energies available 
at flxed target experiments. Here, fragmentation of non- 
strange quarks into kaons is suppressed by an additional 
factor As < 1 , in order to account for the lower probabil- 
ity to generate an ss quark pair rather than a lighter uu 
pair. A value of As = 0.2 was assumed for the HER- 
MES kinematic domain, based on results from deep- 
inelastic muon and neutrino scattering experiments at 
similar center~of-mass energies 85, 86, 87] . This ap- 
proach allows one to compute Pq^^^ {x) directly ac- 
cording to Eq. (|22|l instead of using a Monte Carlo model 
of the fragmentation process. In this analysis, the ele- 
ments of the 2x2 purity matrix were calculated at the 
central kinematic values {x) and (Q^) in each x-hva, us- 
ing the LO parameterizations of the kaon fragmentation 
functions from Ref. [s^ and the CTEQ5L unpolarized 
parton densities from Ref. [60|. 

By computing the purities from parameterizations of 
the fragmentation functions instead of from a Monte 
Carlo model, one does not take into account acceptance 
effects on the measured asymmetries. However, as ex- 
plained in section the effects of the limited polar ac- 
ceptance of the HERMES spectrometer on selected semi- 
inclusive asymmetries are negligible. The asymmetries 
were corrected for azimuthal acceptance effects. The pre- 
cision of the method presented in this section is limited by 
the statistical uncertainties of the measured kaon asym- 
metries, and the knowledge of the strangeness suppres- 
sion factor As. 

Fig. 1221 shows the semi-inclusive K'^ + K~ asymme- 
try on a deuterium target, corrected for effects of QED 
radiation and detector smearing, and azimuthal detector 
acceptance. The same kinematic limits and extraction 
method as described in section were applied. The nu- 
merical values can be found in Tab. IXlTTl in App. B. 
Combined with the inclusive asymmetry Ai ^ shown in 
Fig. the total strange quark polarization [AS'/S'](x) 
was extracted from this data set. The average strange 
quark helicity density 
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FIG. 22: The semi-inclusive Born level asymmetry A^^^^ 
for the total charged kaon flux on a deuterium target, cor- 
rected for instrumental smearing and QED radiative efi'ects. 
The uncertainties are presented analogously to Fig. 1131 



is obtained by multiplication with the unpolarized 
strange quark density s{x) — s{x). The second equality 
in Eq. H45|) holds because the unpolarized strange quark 
density is symmetric by construction for the quark and 
anti-quark in the parameterization employed [20|- The 
resulting average strange quark helicity density is shown 
in Fig. In analogy to the determination of the po- 
larizations of the sea flavors in the previous section, this 
extraction was restricted to the range 0.023 < x < 0.3. 

The uncertainty in the strangeness suppression factor 
contributes the dominant part of the systematic uncer- 
tainty on the extracted values of a;AS'(a;). It was esti- 
mated by varying in the range 0.15 < < 0.3 covered 
by the experimental results [sl, |Q (S^l . This contribu- 
tion to the total systematic uncertainty is shown as the 
unshaded band in Fig. 1231 The remaining contribution 
shown as the shaded band arises from the systematic un- 
certainties on the Born level asymmetries treated in the 
previous section. 

The results from the five flavor fit to the full data set 
on both the proton and deuteron targets are also shown 
in Fig. 1231 A comparision of the first moments of the 
helicity densities in the measured x region provides a 
good measure of the agreement between the two meth- 
ods. It should be noted that these two analyses use the 
same input data. For the five flavor fit, from Tab. IVIlH 
As(5par) = 0.028 ± 0.033 ± 0.009. This is to be com- 
pared with AS'(iso) = 0.129 ± 0.042 ± 0.129 where the 
systematic uncertainty is dominated by that in A^ cited 
above. Considering statistical and systematic uncertain- 
ties, both results are consistent with zero. (After this 
analysis was completed it was learned that the charged 
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FIG. 23: The average strange quark helicity density x ■ 
[As{x) -\- As{x)]/2 from the isoscalar extraction method (full 
points). For comparison the open symbols denote the results 
from a five parameter fit (see text), which are ofi'set horizon- 
tally for presentation. The band in the bottom part gives the 
total systematic uncertainty on the results from the isoscalar 
extraction. The dark shaded area corresponds to the uncer- 
tainty from the input asymmetries, and the open part relates 
to the uncertainty in As. 



kaon multiplicities were consistent with a value somewhat 
larger than this range.) The helicity densities from the 
isoscalar method and from the full five flavor separation 
agree within their uncertainties. The isoscalar method 
yields results that are consistent with a vanishing total 
strange quark helicity density [As(x) -I- As(a;)] = 0. In 
analogy to the previous section, the ^ jndj value for this 
hypothesis is 5.3/7 for the results from the isoscalar ex- 
traction method. Neither the isoscalar nor the five flavor 
measurement provides any indication of a negative po- 
larization of the strange sea. 



VIII. MOMENTS OF HELICITY 
DISTRIBUTIONS 

A. Determination of the Moments 

In the measured region, the n-th moment of the helic- 
ity distributions is given by 



A(")g(Qg) 



E 



Ao 

9 h. 



(46) 

where the quark polarization [Aq/q] (a;,) is assumed to 
be constant in each x-bin i and q(x, Qq) is ob- 

tained from the CTEQ5L LO parameterization of the un- 
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polarized quark densities from Ref. |60j at a fixed scale 
of (3§ = 2.5 GeV^. For simplicity, the first moment 
A(i)g(Qg) is denoted Aq{Ql) hereafter. 

The simultaneous fit of all quark flavors in all x-bins 
yields the correlations between different x-bins as well 
as between different quark flavors, which are taken into 
account for the computation of the uncertainties on the 
moments. Specifically, from the solution of Eq. (|39|l one 
obtains the statistical and systematic covariance matrices 
Vq'^^ and Vq^ for the quark polarizations. The statistical 
and systematic uncertainties on the moments, 6^^^ and 
^Aq ^ obtained from these covariance matrices as 

hi 

X f^^'^dxx''-^ q{x,Ql) f^'^'dxx''-^ q{x,Ql), (47) 

where i and j run over all x-bins. By the inclusion of 
the correlations between different x-bins the total un- 
certainties on the moments are reduced by 22 - 39%, 
depending on the quark flavor. This can be understood 
as a partial compensation of the uncertainty enlargement 
on the asymmetries due to inter-bin event migration in 
the unfolding procedure (cf. section IV Dl and Add. 
The net migration of events from x-bin i into a;~bin j 
causes these bins to be anticorrelated and hence reduces 
the resulting uncertainty in Eq. (|47|l . 

Theoretical predictions are often made in terms of the 
first moments of the isospin singlet, the isovector, and 
the octet combinations, 

AE Au + Au + Ad + AJ+ As + As, (48) 
A93 = Au + Au - (Ad + AJ) , (49) 
A98 = Au + All + Ad + A J - 2( As + As) . (50) 

The resulting values in the measured range 0.023 < a; < 
0.6 are tabulated together with the first and second mo- 
ments of the individual helicity densities in Tab. IVlllI 
For the sea flavors, the measured range is 0.023 < a; < 
0.3. The contributions from 0.3 < a; < 0.6 to these mo- 
ments have been fixed at zero. The valence distributions 
are obtained as Amv = Am — Am and Adv = Ad — Ad. In 
the computation of the uncertainties of these flavor com- 
binations the correlations between the individual helic- 
ity densities are taken into account. The corresponding 
statistical and systematic correlations for the first and 
second moments are listed in Tabs . IXVH and IX VIII in the 
Appendix. 

B. Comments on Extrapolations outside the 
Measured Range 

In order to compute the full moments for the com- 
plete x-range < x < 1, the contributions outside the 
measured region would have to be estimated. For this. 



TABLE VIII: First and second moments of various he- 
licity distributions in the measured range at a scale of 
Ql = 2.5 GeV^ 



Moments in measured range 


Au 


0.601 ± 0.039 ± 0.049 


Au 


-0.002 ± 0.036 ± 0.023 


Ad 


-0.226 ± 0.039 ± 0.050 


Ad 


-0.054 ±0.033 ±0.011 


As 


0.028 ± 0.033 ± 0.009 


All + Au 


0.599 ± 0.022 ± 0.065 


Ad + Ad 


-0.280 ± 0.026 ± 0.057 


AUv 


0.603 ± 0.071 ± 0.040 


Adv 


-0.172 ±0.068 ±0.045 


Au — Ad 


0.048 ± 0.057 ± 0.028 


AS 


0.347 ± 0.024 ± 0.066 


Aqs 


0.880 ± 0.045 ± 0.107 


Aqs 


0.262 ± 0.078 ± 0.045 


A(^'« 


0.142 ±0.009 ±0.011 


A(2)n 


-0.001 ± 0.005 ± 0.002 


A(2)d 


-0.049 ±0.010 ±0.013 


A(2)d 


-0.003 ± 0.004 ± 0.001 


A(2)s 


0.001 ± 0.003 ± 0.001 


A(^'uv 


0.144 ±0.013 ±0.011 


A(2)dv 


-0.047 ±0.012 ±0.012 



one has to rely on models in the kinematic regions where 
no data exist. In particular, the extrapolation towards 
X = is very problematic with many cornpeting models 
making contradictory predictions [Sa . Is^ l90| . This may 
be compared to the situation regarding the behavior of 
the unpolarized structure function -^2(2;, Q^) towards low 
values of x before the HERA data became available (see 
for instance [IJ). 

While in our earlier publication we chose one par- 
ticular model (based on Regge theory) , we do not pursue 
this approach any longer. The reason for this is the very 
strong model dependence which cannot be quantified re- 
liably in terms of an associated systematic uncertainty, as 
some of the models yield diverging integrals in the low-x 
range [9^ . 

To illustrate this more clearly, we performed a st udy 
of the robustness of a global QCD fit similar to 
(LO, scenario 1), which does include preliminary HER- 
MES inclusive data taken on a polarized deuterium tar- 
get 93]. In this study, the fit was repeated with ar- 
tificially offset values for the parameters ry'^ = 0.730 
and 77^ = —0.270 as compared to their default settings 
Vu^ = 0.926 and 7/^^ = —0.341, obtained from the con- 
stants F and D measured in weak neutron and hyperon 
/3~decays. The parameterizations for the helicity densi- 
ties in the QCD fit are constructed such that the param- 
eters r]u^ and rjd^ give the first moments of the helicity 
valence quark distributions, rju.^ = dx Amv(x, Q^) and 
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This modification implies a 21 % change of the value 
of the Bjorken Sum Rule (BJSR) [llil 



h = Aga ^ I dx A(j3(x) 
It) 



9v 



(51) 



from its physical value I3 — \gA/gv\ = 1.267 to = 1.0. 
The unmodified fit agrees well with the obtained helic- 
ity distributions in the measured range as illustrated in 
Fig. 1201 and Tab. IVIII in section IVL^ This consistency 
is also refiected in the good agreement of the first mo- 
ments in the measured range of the valence quark helic- 
ity densities Au^ and Ad^ as well as the isotriplet flavor 
combination Aq^ from the BB LO parameterization with 
the present results from semi-inclusive data, as listed in 
Tab. IIXI It is, however, remarkable that the moments of 
the same quantities computed from the modified fit yield 
almost as good consistency with the measured results, de- 
spite the significant changes in the boundary conditions 
riu^ and 77^^ . 



TABLE IX: Comparison of first moments of the LO QCD fit 
[6^ and a modified version of it (see text) with the measured 
results, in the range 0.023 < a; < 0.6. The values in square 
brackets give the absolute deviation of the moments of the 
QCD fits from the measured SIDIS values in units of their 
combined statistical and systematic uncertainties. 





Original fit 


Modified fit 


Measured SIDIS 




0.601 [0.02 cr] 


0.569 [0.40 cr] 


0.603 ± 0.081 




-0.215 [0.50 cr] 


-0.220 [0.56 cr] 


-0.172 ±0.082 


Aga 


0.816 [0.52 cr] 


0.789 [0.74 cr] 


0.880 ±0.116 



While the QCD fit appears to be stable in the kine- 
matic region where it is rather well constrained by data 
(the data sets entering the fit cover a range 0.005 < 
X < 0.75) it seems to exhibit enough freedom in the un- 
measured regions to allow for a drastic reduction of the 
total integral I3 from 1.267 to 1.0. Table Ixl lists the frac- 
tional contributions to the moments from the original 
and the modified QCD fit in the range < a; < 0.023. In 
order to match the smaller value of /g — 1.0, the low-x 
contributions to all moments are significantly reduced in 
the modified fit. In particular, the low-cc contribution 
from the original fit J^ '^^^dx Aq3{x) — 0.411 reduces to 
J^ '^^^dx Aq'^{x) ~ 0.178 in the modified fit. For com- 
parison, the high-a; extrapolation J^^ dx Aq^{x) — 0.041 

vs. J^^dx Aq';^{x) — 0.033 is only slightly reduced and 
despite the large covered x-range its contribution to the 
total moment is small. 

This study provides information about the "fiexibility" 
of the chosen parameterizations rather than covering the 
full range of competing theoretical models for the extrap- 
olations into the unmeasured kinematic regions. Yet, it 
demonstrates the arbitrariness of these extrapolations, in 
particular to the small x range. In the measured region, 
our experimental result on Aq^ agrees within 0.52 a with 



TABLE X: Fractional contributions from the low-a; range 
(0 < X < 0.023) to the first moments obtained from the orig- 
inal and modified global QCD fit (see text). 





Original fit 


Modified fit 




31% 


18% 


Adv 


36% 


18% 


Aga 


32% 


18% 



the integral of this quantity from the original QCD fit 
over the same x-range (see Tab. IIX|I . We nevertheless 
refrain from interpreting this as an experimental confir- 
mation of the BJSR Eq. (|51|l . as the above example has 
shown that one can obtain a similar level of agreement 
with a 21 % reduction of the predicted value, thereby 
rendering a sound estimate of the associated systematic 
uncertainty impossible. For the same reason we restrict 
all results for the first and second moments to the experi- 
mentally covered a;-region, albeit the weight of the low-x 
extrapolation becomes smaller for the higher moments. 

We finally remark that in j96l l a reverse approach to 
this problem is presented, where the authors assume the 
validity of the BJSR Eq. and use experimental data 
at X > 0.005 to restrict the exponent A = —0.40 ±0.31 in 
an assumed power-like behaviour of the polarized struc- 
ture function gi{x) ~ x^, based on Regge theory. 



C. Comparison with other Results 

SMC is the only other experiment that has published 
results from SIDIS on the quark helicity distributions in 
the nucleon 48]. Due to limited statistics and the lack 
of discrimination between different hadron types, SMC 
could extract quark helicity distributions only under the 
assumption of SU{3) flavor symmetry for the sea quark 
flavors, i.e. Au{x) = Ad{x) = As{x) = As{x). It is only 
because of this assumption that their precision in sea 
quark polarization is comparable to that of the present 
work. 

The first and second moments for the valence and u 
sea quark helicity densities in the measured range of the 
HERMES data are compared in Table [XTI to the results 
from the SMC experiment, which were integrated over 
the same kinematic region. It should be noted that SMC 
evaluated their moments according to a somewhat sim- 
pler expression 97] 



A(")g(g^) - ^(6+1 - e.) Aq(x., Q; 



(52) 



where the notations have been adapted to match the def- 
inition in Eq. (|46l) . The integration procedure assumes 
constant helicity densities Aq{x) over the entire width of 
an x-bin. In particular towards the upper experimental 
limit on x, where the x-bins are widest, the values de- 
termined for the helicity densities Auv(x) and Adv{x) in 
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TABLE XI: Comparison of first and second moments in the 
measured range from this analysis with results from the SMC 
experiment. The SMC values were extrapolated to the same 
value of Ql = 2.5 GeV^ and integrated over the HERMES 



a;-range. 





HERMES 


SMC 




0.603 ± 0.071 ± 0.040 


0.614 ± 0.082 ± 0.068 


Adv 


-0.172 ±0.068 ±0.045 


-0.334 ±0.112 ±0.089 


Am 


-0.002 ± 0.036 ± 0.023 


0.015 ±0.034 ±0.024 


A(2)«, 


0.144 ±0.013 ±0.011 


0.152 ±0.016 ±0.013 


A(2)d, 


-0.047 ±0.012 ±0.012 


-0.056 ±0.026 ±0.015 



some cases violate the positivity limits given by the un- 
polarized parton densities at the upper x-bin limits (see 
Fig. 2 in j^). This causes larger values for the moments 
when computed according to Eq. (|52|) as compared to 
Eq. (|46l) . Nevertheless the results from the two experi- 
ments are in good agreement within their combined un- 
certainties, while the present data set has an improved 
precision. When comparing the precision of the results 
from both experiments, note that the HERMES results 
are free of symmetry assumptions regarding the sea quark 
helicity densities and that they are based on a new un- 
folding technique to account for the proper propagation 
of uncertainties while handling QED radiative and in- 
strumental smearing effects. 



IX. SUMMARY 

This paper describes the most precise available mea- 
surements of semi-inclusive asymmetries in polarized 
deep inelastic electron/positron scattering. Results from 
both inclusive and semi-inclusive measurements are pre- 
sented. The unique semi-inclusive data are particularly 
important for the unbiased extraction of sea quark helic- 
ity densities in the nucleon. 

Data were collected on longitudinally polarized atomic 
hydrogen and deuterium gas targets. Good particle iden- 
tification in the HERMES spectrometer allows hadrons 
that are coincident with the scattered lepton to be sepa- 
rated into pion and kaon samples. This gives sensitivity 
to the individual quark and antiquark helicity distribu- 
tion functions for light and strange quark flavors. 

An unfolding technique — new to this type of measure- 
ment — was applied for instrumental smearing and QED 
radiative effects. It takes into account event migrations 
between bins. This algorithm provides more rigorous, 
yet larger, estimates of the inflation of the uncertainties 
on the Born asymmetries with respect to the measured 
asymmetries than correction methods previously applied. 
The unfolding procedure also yields previously unavail- 
able estimates of the statistical correlations between dif- 
ferent kinematic bins. All asymmetries were found to be 
positive and increase with x, with the exception of the 



K~ asymmetry which is consistent with zero over the 
entire measured range. The present data on undiscrim- 
inated hadrons agree well with earlier measurements by 
SMC at higher Q^, albeit with much improved precision. 

A "leading order" QCD analysis that relies on the tech- 
nique of flavor tagging was used to extract helicity dis- 
tribution functions. The results presented here allow for 
the first time the independent determination of five out 
of six quark polarizations in the nucleon. Quark po- 
larizations [lS.q/q][x) were obtained using purities that 
were calculated in a Monte Carlo based on parameteri- 
zations of unpolarized parton densities and a modeling 
of hadron multiplicities measured at HERMES. The po- 
larization [lS.u/u]{x) was found to be positive and rising 
over the entire range in a;, while [lS.d/(I\{x) is negative. 
These first results on the individual sea quark polariza- 
tions [/!^u/u]{x), [AJ/J](x), and [As/s](a;) are consistent 
with zero. 

Furthermore, the "LO" approach presented here re- 
sults in helicity quark distributions that are fully consis- 
tent with global QCD fits to world data on inclusive DIS 
asymmetries. While the BB fit favors a slightly negative 
strange quark helicity density As(a;), the data yield a 
small positive result. However, it should be noted that 
those particular fits invoke SU (3) symmetry to relate the 
triplet as and octet ag axial couplings to the weak decay 
constants F and D, assume a SU (3) symmetric sea, and 
require a model for the gluon helicity density l^g{x). Be- 
cause of these assumptions and taking into account the 
combined statistical and systematic uncertainties on the 
present experimental result, no significant discrepancy 
remains. However, the results from the present semi- 
inclusive scattering data might reveal the possible biases 
of global fits for certain quark flavors. Within the exper- 
imental uncertainty all obtained quark helicity densities 
Ag(x) are in good agreement with the most recent pa- 
rameterizations. 

To confirm the results for the helicity distribution of 
the strange SG0(, Oil different isoscalar analysis leading to 
As(a;) + As(a;) was performed. This technique does 
not depend on a Monte Carlo modehng of the strange 
fragmentation functions. Instead, parameterizations of 
strange and non-strange fragmentation functions mea- 
sured at e^e~ colliders at higher energies were utilized. 
The results for the two analyses are in reasonable agree- 
ment. 

The first direct experimental extraction of the helic- 
ity density asymmetry Am(x) — /S.d{x) in the light quark 
sea, which is predicted to be non zero by many models in 
analogy to the unpolarized sector {u{x) ^ d{x)), does not 
establish broken SU{2) flavor symmetry. The data disfa- 
vor the substantial positive asymmetry predicted by the 
xQSM model and are consistent with the small negative 
asymmetry characteristic of the meson cloud model. 

Moments were computed in the measured kinematic 
range for the quark helicity distributions and singlet and 
non-singlet flavor combinations. The moments of the 
valence and sea quark helicity densities Auv, Adv, and 
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Au agree with previous results from the SMC experi- 
ment. In the measured range, the non-singlet flavor 
combination Aq^ is in good agreement with the same 
quantity computed from global QCD fits, which in turn 
by construction fulfill the Bjorken sum rule. However, 
because extrapolations into the unmeasured kinematic 
regions are not sufficiently constrained by the data re- 
ported here, an experimental confirmation of this funda- 
mental sum rule was not possible. For the same reason 
we do not quote a total value for the singlet quantity 
AE, which gives the total contribution from quark spins 
to the spin of the nucleon. In the measured range, a value 
of AS = 0.347 ± 0.024 ± 0.066 was obtained. 

In conclusion, HERMES has made detailed semi-inclu- 
sive measurements of polarization asymmetries in deep 
inelastic lepton scattering. These measurements allow 
quark helicity distribution functions to be extracted with 
fewer and different model assumptions than in previous 
inclusive measurements. 
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APPENDIX A: UNFOLDING OF RADIATIVE 
AND DETECTOR SMEARING EFFECTS 

Corrections to the asymmetries for higher order QED 
and detector smearing effects were carried out using an 
unfolding algorithm as indicated in section IV DI How- 
ever, the algorithm is complex, because the known ab- 
solute cross sections of the spin-dependent background 
processes must be normalized to the data by compar- 
ing simulated and measured unpolarized yields based 
on world data for ^2- The procedure uses two sets 
of Monte Carlo data. Born data were generated with 
the GRSV2000 (LO, standard scenario) spin-dependent 
parton distributions and the LUND fragmentation 
model implemented in Jetset 7.4 Within the ac- 

ceptance of the HERMES spectrometer < 170mrad, 
40mrad < \9y\ < 140mrad) an equivalent number of 
11.5 M DIS events were generated on hydrogen and deu- 
terium targets. The second Monte Carlo data set in- 
cluded, in addition, internal radiative effects and a sim- 



ulation of the HERMES spectrometer. Radiative effects 
were simulated with RADGEN and the spectrome- 
ter simulation was done using GEANT [§3 ■ An equiva- 
lent number of 5.4 M DIS events was available for analysis 
for each target. The total MC statistics are sufficiently 
large that the uncertainty on the Born asymmetries from 
this source is smah (cf. Tabs. IXIll IXHll and IXIVjl . 

Based on the latter Monte Carlo data, matrices 
n^^^j(i,j) with dimensions nx x (ns + 1) were con- 
structed for the parallel (=^) and anti-parallel (^) spin- 
states, which describe the count rates that fall in both 
bin j of Born kinematics and bin i of e2Cperimental kine- 
matics. The experimental kinematics reflect the radia- 
tive and instrumental effects. The bins j = 1, . . . , ub 
describe the migration of DIS events, where both the 
experimental and the Born kinematics are within the ac- 
ceptance. Here, acceptance refers to both the geometri- 
cal acceptance of the spectrometer and the phase space 
defined by the DIS and SIDIS cuts. The additional bin 
j = contains background rates that feed into the exper- 
imental bins through QED radiative and detector smear- 
ing effects. For example, in the case of the inclusive data 
sample on the proton, the background stems from elastic 
and inelastic scattering events that are radiated (through 
QED effects) or smeared (by interactions in the detector) 
into the acceptance. The experimental rates in the par- 
allel and anti-parallel spin states in bin i are given by the 
sum (-)(«) = J2]=o 'n'^(^)ihj)- The corresponding 



Born rates, "■!^(-^j(i)i iii each spin state and bin j were 
calculated from the Born Monte Carlo data. 

These data, normalized with respect to each other pro- 
vide access to the migration matrices '5^^. — ^(z, j), which 
are given by 

'^-(-)(*'^)-9a|^^^(jr 

1 = 1,..., nx, j = 1, . . . , n_B . (Al) 

In terms of n,^^^^ (i, j) and "^^^^j (j)j these matrices can 
be written as 



(A2) 



provided the first derivatives are constant. The matrices 
describe the kinematical migration inside of the 
acceptance. Backgrounds are taken into account through 
the normalization with the Born rates n?, , ^ and as will 

be seen below in the j = row of «^(^)(*jj)- The 
^^{^) matrices are insensitive to the Monte Carlo model 
of the Born distributions, because both the numerator 
and the denominator scale with the relative number of 
Born events generated in bin number j- 

The unknown Born rates B^^^^{j) are related to the 
measured asymmetries (Eq. ()30|l corrected for the az- 
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imuthal acceptance effects) through the migration matri- 
ces, the mipolarized experimental rates {i) = + 

n^(i) and the spin-dependent background np{i,0) = 

Af{i)nu{i)~np{i,0), i = 1, . . . , . (A3) 

Here the unpolarized Born rate B^{j) + B^{j) = Bu{j) 
is known from previous experiments and incorporated in 
the Monte Carlo simulation: Bu{j) = n^{j) = +n^. 
The background term in the sum was moved to the right 
hand side of the equation np{i,0) — S^{i,0)B^{0) — 
S~{i,Q)B^{0). The rate B^{j) may now be eliminated 
in favor of B^{j) and Bu{j). 



-np{i,0) + ^S^{i,j)n^{j), i = 1, . . . , nx ■ 

i=i 



(A4) 



The Born asymmetry is found by solving Eq. (|A4|) for 
B^{j) and substituting into 



, 2B^ij)-B^ 
"^11 - B^ij) 



(A5) 



The final expression for the unfolded asymmetry is then: 



Af{i) (0 - np{i, 0) + £ S^{i, k) n^{k) 



k=l 



(A6) 
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FIG. 24: Matrices nu{i,j) = n^{i, j)+n^{i, j) for DIS events 
and SIDIS vr"*" events on the proton. The binning shown cor- 
responds to the 9 bins in x that were used in the asymmetry 
and Aq analysis (Tab. [3 . See the text for details. 



for j = 1, . . . , ns, where S'{i,j) is the square submatrix 
without the column j = 0of5(j,j) = S^{i,j)+S^{i,j). 
General ly th e inverse of the migration matrix may be ill- 
defined JIOQI . However, in this case the migration due to 
QED and detector effects is sufficiently small that the 
matrix inversion yields satisfactory results. 

Fig- ISl presents the matrices nu{i,j) for DIS events 
and SIDIS tt^ events on the proton calculated in the x- 
Bjorken bins of the analysis presented in this paper. The 
background that migrates into the acceptance (j = 0) 
is seen to be large in case of the DIS events. Only 
little background is present in the case of the SIDIS 
events, where the hadron-tag rejects elastic events that 
are present in the inclusive case. The inter-bin migration 
is of similar size in both data-samples. QED radiative ef- 
fects, specifically initial and final-state bremsstrahlung, 
cause migration to smaller values of observed x only. 
Multiple scattering and finite resolution effects in the de- 
tector also increase but may in some cases decrease the 
observed x with respect to the Born x. 

The statistical covariance matrix of the Born asymme- 
try follows from Eq. (jA6|l 



ViAf{j),Af{k))^ 



(A7) 



ii=l 12 = 1 



where 2?(j, i) is the dilution matrix. 



2?(j,0 



dAfjj) 
dAf{i) 



(A8) 



These expressions unambiguously determine the statisti- 
cal uncertainties on the Born asymmetries, unlike pre- 
vious methods that compute the corrections bin-by- 
bin. In a similar manner, systematic uncertainties due 
to beam and target polarization, as well as for the az- 
imuthal acceptance correction were determined for the 
measured aysmetries. The systematic covariances from 
these sources follow via the dilution matrix. For the beam 
and target polarization measurements the uncertainty in 
the Born asymmetries is 

Vp,,,,(Af(j),Af(fc)) = 

nx nx 
ii = l 12 = 1 

(A9) 

The double sum over the experimental uncertainties is 
a consequence of the complete correlation of the beam 
and target polarization measurements among the x bins. 
In the case of the azimuthal acceptance correction, the 
systematic uncertainty was assumed to be uncorrelated 
between the bins. 



APPENDIX B: TABLES OF RESULTS 
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TABLE XII: Inclusive and semi-inclusive Born level asymmetries on the proton target. 



{x) (Q^)/GeV^ Ai,p ±stat. ±syst. ±MC 



0.033 


1.22 


0.0996 


0.0176 


0.0091 


0.0037 


0.048 


1.45 


0.1102 


0.0175 


0.0093 


0.0043 


0.065 


1.68 


0.1131 


0.0181 


0.0088 


0.0050 


0.087 


1.93 


0.1941 


0.0208 


0.0143 


0.0063 


0.118 


2.34 


0.2366 


0.0205 


0.0178 


0.0061 


0.166 


3.16 


0.2819 


0.0216 


0.0208 


0.0054 


0.240 


4.54 


0.3854 


0.0239 


0.0283 


0.0045 


0.340 


6.56 


0.4760 


0.0430 


0.0344 


0.0054 


0.447 


9.18 


0.6102 


0.0750 


0.0467 


0.0059 



(x) (Q^)/GeV^ (z) aC istat. ±syst. ±MC A^p ±stat. ±syst. ±MC 



0.034 


1.21 


0.356 


0.1097 


0.0299 


0.0077 


0.0067 


0.0724 


0.0341 


0.0054 


0.0073 


0.048 


1.44 


0.375 


0.1640 


0.0316 


0.0113 


0.0075 


0.1262 


0.0365 


0.0087 


0.0085 


0.065 


1.72 


0.386 


0.1361 


0.0329 


0.0096 


0.0083 


0.0901 


0.0382 


0.0057 


0.0094 


0.087 


2.06 


0.393 


0.2075 


0.0370 


0.0147 


0.0095 


0.0953 


0.0442 


0.0066 


0.0111 


0.118 


2.58 


0.395 


0.3011 


0.0356 


0.0212 


0.0088 


0.1860 


0.0445 


0.0135 


0.0109 


0.166 


3.52 


0.390 


0.2851 


0.0384 


0.0209 


0.0082 


0.1770 


0.0498 


0.0127 


0.0104 


0.239 


5.03 


0.388 


0.4223 


0.0456 


0.0292 


0.0075 


0.2360 


0.0639 


0.0165 


0.0106 


0.338 


7.09 


0.377 


0.4046 


0.0914 


0.0291 


0.0101 


0.5696 


0.1339 


0.0387 


0.0157 


0.448 


9.76 


0.364 


0.7586 


0.1756 


0.0519 


0.0122 


0.4957 


0.2715 


0.0341 


0.0231 


(x) 


(Q2) / GcV^ 






±stat. 


zbsyst. 


±MC 




itstat. 


zbsyst. 


±MC 


0.033 


1.22 


0.364 


0.0800 


0.0353 


0.0058 


0.0077 


0.0675 


0.0388 


0.0053 


0.0081 


0.047 


1.50 


0.416 


0.1336 


0.0387 


0.0091 


0.0088 


0.1450 


0.0431 


0.0095 


0.0092 


0.064 


1.87 


0.449 


0.0829 


0.0408 


0.0067 


0.0093 


0.0649 


0.0461 


0.0039 


0.0100 


0.087 


2.38 


0.471 


0.2312 


0.0459 


0.0157 


0.0102 


0.0714 


0.0536 


0.0047 


0.0109 


0.118 


3.08 


0.487 


0.3163 


0.0458 


0.0212 


0.0093 


0.0754 


0.0547 


0.0057 


0.0101 


0.166 


4.22 


0.490 


0.3017 


0.0525 


0.0201 


0.0090 


0.1572 


0.0645 


0.0105 


0.0102 


0.238 


5.83 


0.504 


0.2784 


0.0695 


0.0197 


0.0092 


0.2696 


0.0889 


0.0187 


0.0108 


0.337 


7.97 


0.506 


0.5566 


0.1530 


0.0373 


0.0138 


0.3461 


0.1995 


0.0233 


0.0177 


0.449 


10.49 


0.496 


0.8651 


0.3185 


0.0558 


0.0175 


0.4490 


0.4343 


0.0352 


0.0270 



[1] M. Alguard ct al. (E80), Pliys. Rev. Lett. 37, 1261 
(1976). 

[2] G. Baum et al. (E130), Phys. Rev. Lett. 51, 1135 (1983). 
[3] R. L. Jaffe and A. Manohar, Nucl. Phys. B 337, 509 

(1990). 

[4] K. Suzuki and W. Wcise, Nucl. Phys. A 634, 141 (1998). 
[5] J. Ashman ot al. (EMC), Phys. Lett. B 206, 364 (1988). 
[6] J. Ashman et al. (EMC), Nucl. Phys. B 328, 1 (1989). 
[7] G. Altarelh, R. D. Ball, F. Forte, and G. Ridolfi, 

Nucl. Phys. B 496, 337 (1997). 
[8] R. Mertig and W. L. van Neerven, Z. Phys. C 70, 637 

(1996). 

[9] W. Vogelsang, Phys. Rev. D 54, 2023 (1996). 

[10] B. Adeva et al. (SMC), Phys. Rev. D 58, 112002 (1998). 

[11] F. Halzen and A. Martin, Quarks and Leptons: An In- 
troductory Course in Modem Particle Physics (John 
Wiley and Sons, 1984). 

[12] R. Roberts, The Structure of the Proton: Deep Inelastic 
Scattering (Cambridge University Press, 1990). 

[13] A. Thomas and W. Weise, The Structure of the Nucleon 
(Wiley- VCH, 2001). 



[14] C. Callan and D. Gross, Phys. Rev. Lett. 22, 156 (1969). 
[15] K. Abe et al. (E143), Phys. Rev. Lett. 76, 587 (1996). 
[16] K. Abe et al. (E154), Phys. Lett. B 404, 377 (1997). 
[17] P. Anthony et al. (E155), Phys. Lett. B 553, 18 (2003). 
[18] M. Gliick, E. Reya, M. Stratmann, and W. Vogelsang, 

Phys. Rev. D 63, 094005 (2001). 
[19] J. Buon and K. Steffen, Nucl. Inst. & Meth. A 245, 248 

(1986). 

[20] D. Barber et al., Nucl. Inst. & Meth. A 329, 79 (1993). 
[21] D. Barber et al., Nucl. Inst. & Meth. A 338, 166 (1994). 
[22] M. Beckmann et al., Nucl. Inst. & Meth. A 479, 334 

(2002). 

[23] A. Sokolov and I. Ternov, Sov. Phys. Doklady 8, 1203 
(1964). 

[24] M. Diiren, Habilitation Thesis, Friedrich-Alexander- 
Universitat Erlangen-Niirnberg (1995), DESY HER- 
MES-95-02. 

[25] D. DeSchepper et al., Nucl. Inst. & Meth. A 419, 16 

(1998). 

[26] A. Nass et al., Nucl. Inst. & Meth. A 505, 633 (2003). 
[27] C. Baumgarten et al., Nucl. Inst. & Meth. A 482, 606 



33 



TABLE XIII: Inclusive and semi-inclusive Born level asymmetries on the deuterium target. 



{X) 


w / / 'jev 








isy st . 












0.033 


1.22 




0.0203 


0.0078 


0.0015 


0.0034 










0.048 


1.45 




0.0248 


0.0080 


0.0017 


0.0040 










0.065 


1.69 




0.0396 


0.0085 


0.0023 


0.0048 










0.087 


1.95 




0.0440 


0.0100 


0.0031 


0.0061 










0.118 


2.35 




0.0777 


0.0099 


0.0054 


0.0060 










0.166 


3.18 




0.1137 


0.0107 


0.0081 


0.0054 










0.240 


4.55 




0.1621 


0.0121 


0.0114 


0.0046 










0.339 


6.58 




0.2932 


0.0228 


0.0195 


0.0057 










0.446 


9.16 




0.3161 


0.0412 


0.0236 


0.0065 










/t\ 
\X) 


W / / ^ 




Ah + 




itsystr. 




Ah 




ztsystj . 


ZIZiVl\_' 


0.033 


1.21 


0.355 


0.0080 


0.0146 


0.0007 


0.0068 


-0.0125 


0.0162 


0.0013 


0.0071 


0.048 


1.44 


0.374 


0.0112 


0.0156 


0.0018 


0.0077 


0.0074 


0.0174 


0.0014 


0.0084 


0.065 


1.73 


0.384 


0.0484 


0.0162 


0.0028 


0.0085 


0.0380 


0.0185 


0.0022 


0.0095 


0.087 


2.07 


0.391 


0.0754 


0.0185 


0.0043 


0.0101 


0.0179 


0.0212 


0.0033 


0.0113 


0.118 


2.60 


0.394 


0.0350 


0.0179 


0.0038 


0.0094 


0.0739 


0.0213 


0.0040 


0.0108 


0.166 


3.56 


0.392 


0.1326 


0.0194 


0.0087 


0.0087 


0.0775 


0.0245 


0.0065 


0.0106 


0.238 


5.04 


0.388 


0.1469 


0.0237 


0.0104 


0.0081 


0.1712 


0.0315 


0.0103 


0.0107 


0.338 


7.12 


0.382 


0.2372 


0.0504 


0.0151 


0.0115 


0.3001 


0.0700 


0.0175 


0.0168 


0.446 


9.61 


0.380 


0.1901 


0.0995 


0.0149 


0.0147 


0.1499 


0.1481 


0.0128 


0.0242 


/-r\ 


\v / / '-JC V 


\2/ 




-l-e-t-at 
ZLObcLL. 


isySL. 








zbsystj . 


ZlZiVlVy 


0.033 


1.22 


0.353 


-0.0172 


0.0175 


0.0011 


0.0076 


-0.0113 


0.0183 


0.0014 


0.0078 


0.047 


1.50 


0.405 


0.0180 


0.0192 


0.0022 


r\ f\f\c>r\ 

0.0089 


—0.0231 


0.0203 


0.0012 


0.0093 


0.064 


1.87 


0.437 


0.0130 


0.0201 


0.0016 


0.0095 


0.0457 


0.0218 


0.0028 


0.0100 


0.087 


2.36 


0.458 


0.0449 


0.0226 


0.0029 


0.0106 


0.0056 


0.0245 


0.0020 


0.0110 


0.118 


3.07 


0.472 


0.0966 


0.0223 


0.0061 


0.0096 


0.0884 


0.0249 


0.0045 


0.0103 


0.165 


4.18 


0.479 


0.1207 


0.0257 


0.0079 


0.0093 


0.0144 


0.0298 


0.0035 


0.0105 


0.238 


5.80 


0.488 


0.1089 


0.0343 


0.0073 


0.0097 


0.2039 


0.0413 


0.0116 


0.0111 


0.338 


7.93 


0.494 


0.3179 


0.0815 


0.0202 


0.0153 


0.3860 


0.0988 


0.0209 


0.0186 


0.446 


10.24 


0.503 


0.0856 


0.1695 


0.0115 


0.0215 


-0.1323 


0.2159 


0.0182 


0.0282 



{x) (Q^) / GeV^ {z) Af^i ±stat. ±syst. ±MC ±stat. ±syst. ±MC 



0.033 


1.22 


0.383 


0.0048 


0.0479 


0.0022 


0.0202 


-0.0471 


0.0597 


0.0039 


0.0212 


0.048 


1.50 


0.424 


0.0171 


0.0496 


0.0043 


0.0232 


0.0312 


0.0661 


0.0041 


0.0246 


0.065 


1.86 


0.457 


0.1469 


0.0504 


0.0083 


0.0245 


0.0097 


0.0701 


0.0051 


0.0262 


0.086 


2.33 


0.484 


0.1220 


0.0561 


0.0079 


0.0270 


-0.0554 


0.0811 


0.0046 


0.0290 


0.118 


3.08 


0.489 


0.0399 


0.0534 


0.0046 


0.0239 


0.0292 


0.0830 


0.0029 


0.0274 


0.165 


4.23 


0.493 


0.1436 


0.0593 


0.0104 


0.0229 


0.0722 


0.0993 


0.0069 


0.0279 


0.238 


5.81 


0.503 


0.1445 


0.0773 


0.0120 


0.0235 


0.0871 


0.1411 


0.0067 


0.0313 


0.336 


7.76 


0.516 


0.4389 


0.1747 


0.0270 


0.0373 


-0.2504 


0.3422 


0.0202 


0.0621 


0.448 


10.20 


0.510 


0.4641 


0.3692 


0.0400 


0.0536 


1.4585 


0.7001 


0.0859 


0.1010 



(x) (Q2)/GeV^ {z) Af^^+'' ±stat. ±syst. ±MC 



0.033 


1.22 


0.372 


-0.0161 


0.0372 


0.0016 


0.0148 




0.047 


1.50 


0.415 


0.0257 


0.0395 


0.0036 


0.0168 




0.064 


1.87 


0.446 


0.1064 


0.0408 


0.0049 


0.0182 




0.086 


2.36 


0.470 


0.0709 


0.0460 


0.0042 


0.0201 




0.118 


3.11 


0.475 


0.0395 


0.0447 


0.0030 


0.0180 




0.165 


4.26 


0.481 


0.1226 


0.0506 


0.0095 


0.0177 




0.238 


5.87 


0.489 


0.1347 


0.0671 


0.0103 


0.0189 




0.336 


7.81 


0.504 


0.2945 


0.1540 


0.0183 


0.0316 




0. 1 17 


il).2() 


O.IOO 


I). 6517 


l).:!251 


0.0 191 


l).()l()7 
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TABLE XIV: The quaxk polarizations [IS.q/q\{x), and the quark helicity densities x ■ Ag(a;, Qq) evolved to Qo = 2.5 GeV^. The 
systematic uncertainty due to the purities ([±Pttr]) is included in the total systematic uncertainty (±sys). 



{x) 


Am/w 


±stat 


±sys 


[±Pur] 


±MC 


X ■ Am 


±stat 


±sys 


[±Pur] 


±MC 


0.033 


0.0855 


0.1180 


0.0856 


0.0837 


0.0129 


0.0296 


0.0408 


0.0300 


0.0290 


0.0045 


0.048 


0.1368 


0.1107 


0.0217 


0.0037 


0.0077 


0.0515 


0.0417 


0.0101 


0.0014 


0.0029 


0.065 


0.1913 


0.0978 


0.0356 


0.0264 


0.0051 


0.0785 


0.0401 


0.0160 


0.0108 


0.0021 


0.087 


0.4864 


0.0909 


0.0784 


0.0703 


0.0057 


0.2185 


0.0408 


0.0357 


0.0316 


0.0026 


0.118 


0.5086 


0.0774 


0.0633 


0.0501 


0.0057 


0.2525 


0.0385 


0.0317 


0.0249 


0.0029 


0.166 


0.4731 


0.0757 


0.0436 


0.0267 


0.0045 


0.2623 


0.0420 


0.0244 


0.0148 


0.0025 


0.239 


0.4445 


0.0855 


0.0364 


0.0109 


0.0028 


0.2652 


0.0510 


0.0217 


0.0065 


0.0017 


0.339 


0.5805 


0.0650 


0.0558 


0.0050 


0.0091 


0.3241 


0.0363 


0.0313 


0.0028 


0.0051 


0.447 


0.7272 


0.1087 


0.0684 


0.0139 


0.0084 


0.3121 


0.0467 


0.0294 


0.0060 


0.0036 


{x) 


Ad/d 


±stat 


±sys 


[±Pur] 


±MC 


X ■ Ad 


istat 


±sys 


[±Pur] 


±MC 


0.033 


-0.1236 


0.1529 


0.0390 


0.0199 


0.0153 


-0.0337 


0.0416 


0.0128 


0.0054 


0.0042 


0.048 


0.0588 


0.1452 


0.0543 


0.0435 


0.0115 


0.0167 


0.0411 


0.0166 


0.0123 


0.0032 


0.065 


-0.1336 


0.1307 


0.0362 


0.0189 


0.0122 


-0.0394 


0.0385 


0.0120 


0.0056 


0.0036 


0.087 


-0.2572 


0.1303 


0.0495 


0.0165 


0.0190 


-0.0789 


0.0400 


0.0159 


0.0051 


0.0058 


0.118 


-0.4876 


0.1185 


0.0841 


0.0592 


0.0210 


-0.1552 


0.0377 


0.0271 


0.0188 


0.0067 


0.166 


-0.0918 


0.1337 


0.0675 


0.0091 


0.0169 


-0.0296 


0.0431 


0.0219 


0.0029 


0.0055 


0.239 


-0.5218 


0.1646 


0.0822 


0.0015 


0.0125 


-0.1536 


0.0484 


0.0242 


0.0004 


0.0037 


0.339 


-0.2799 


0.1988 


0.1694 


0.0166 


0.0400 


-0.0628 


0.0446 


0.0382 


0.0037 


0.0090 


0.447 


—0.8133 


0.4074 


0.2454 


0.0564 


0.0404 


—0.1158 


0.0580 


0.0349 


0.0080 


0.0057 


(a;) 


Au/u 


±stat 


±sys 


[±Pur] 


±MC 


X ■ Am 


istat 


±sys 


[±Pur] 


±MC 


0.033 


0.3382 


0.3342 


0.2189 


0.2169 


0.0520 


0.0437 


0.0432 


0.0283 


0.0280 


0.0067 


0.048 


0.2484 


0.3677 


0.0471 


0.0074 


0.0448 


0.0288 


0.0426 


0.0058 


0.0009 


0.0052 


0.065 


0.0166 


0.3938 


0.1200 


0.1132 


0.0512 


0.0017 


0.0403 


0.0126 


0.0116 


0.0052 


0.087 


-0.7151 


0.4585 


0.3508 


0.3468 


0.0546 


-0.0624 


0.0400 


0.0307 


0.0302 


0.0048 


0.118 


-0.8989 


0.5391 


0.3395 


0.3324 


0.0626 


-0.0621 


0.0372 


0.0235 


0.0229 


0.0043 


0.166 


-0.9022 


0.8403 


0.3491 


0.3307 


0.0815 


-0.0432 


0.0402 


0.0168 


0.0158 


0.0039 


0.239 


1.4742 


1.6410 


0.3868 


0.2844 


0.1302 


0.0446 


0.0496 


0.0117 


0.0086 


0.0039 


{x) 


Ad/d 


±stat 


±sys 


[±Pur] 


±MC 


X ■ Ad 


±stat 


±sys 


[±Pur] 


±MC 


0.033 


-0.2281 


0.2819 


0.0380 


0.0198 


0.0200 


-0.0338 


0.0417 


0.0067 


0.0029 


0.0030 


0.048 


-0.6238 


0.2921 


0.1076 


0.1033 


0.0387 


-0.0862 


0.0404 


0.0151 


0.0143 


0.0053 


0.065 


0.0174 


0.2847 


0.0513 


0.0457 


0.0446 


0.0022 


0.0367 


0.0068 


0.0059 


0.0058 


0.087 


-0.2239 


0.3103 


0.0605 


0.0550 


0.0482 


-0.0267 


0.0370 


0.0073 


0.0065 


0.0057 


0.118 


0.5412 


0.3144 


0.1621 


0.1588 


0.0449 


0.0577 


0.0335 


0.0173 


0.0169 


0.0048 


0.166 


-0.9546 


0.4561 


0.0734 


0.0545 


0.0590 


-0.0828 


0.0396 


0.0064 


0.0047 


0.0051 


0.239 


0.4523 


0.8380 


0.1470 


0.0323 


0.1075 


0.0237 


0.0439 


0.0077 


0.0017 


0.0056 


{x) 


As/s 


zbstat 


±sys 


[±Pur] 


±MC 


X ■ As 


itstat 


±sys 


[±Pur] 


±MC 


0.033 


0.4734 


0.7492 


0.1871 


0.1659 


0.0596 


0.0317 


0.0502 


0.0134 


0.0111 


0.0040 


0.048 


0.6071 


0.6361 


0.1952 


0.1908 


0.0462 


0.0365 


0.0382 


0.0118 


0.0115 


0.0028 


0.065 


-0.0537 


0.5805 


0.0441 


0.0329 


0.0550 


-0.0029 


0.0312 


0.0024 


0.0018 


0.0030 


0.087 


-0.1243 


0.6248 


0.1422 


0.1366 


0.0661 


-0.0059 


0.0295 


0.0067 


0.0065 


0.0031 


0.118 


-0.3359 


0.6516 


0.0756 


0.0689 


0.0772 


-0.0133 


0.0258 


0.0030 


0.0027 


0.0031 


0.166 


1.3956 


0.8851 


0.0912 


0.0252 


0.0780 


0.0418 


0.0265 


0.0027 


0.0008 


0.0023 


0.239 


-1.2674 


1.5039 


0.3911 


0.3781 


0.1220 


-0.0230 


0.0273 


0.0071 


0.0069 


0.0022 
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TABLE XV: The light sea flavor heUcity asymmetry evolved to Qq = 2.5 GeV^. The systematic uncertainty due to the purities 
([iPwr]) is included in the total systematic uncertainty (±sys). 



{x} 


X ■ (Aw - Ad) 


±stat 


±sys 


[±Pur] 


±MC 


0.033 


0.0748 


0.0653 


0.0252 


0.0242 


0.0093 


0.048 


0.1133 


0.0670 


0.0153 


0.0132 


0.0102 


0.065 


-0.0005 


0.0639 


0.0180 


0.0174 


0.0106 


0.087 


-0.0355 


0.0655 


0.0370 


0.0367 


0.0100 


0.118 


-0.1198 


0.0606 


0.0404 


0.0399 


0.0088 


0.166 


0.0397 


0.0709 


0.0221 


0.0211 


0.0088 


0.239 


0.0241 


0.0885 


0.0178 


0.0077 


0.0097 



TABLE XVI: The statistical (pstat) and systematic (psyst) correlations of the first moments of the extracted helicity distributions 
in the measured x-range. 



Pstat 


Am 


Ad 


Am 


Ad 


As 


psyst 


Am 


Ad 


Am 


Ad 


As 


Am 


1.000 


-0.564 


-0.827 


0.346 


0.075 


Am 


1.000 


-0.732 


0.577 


-0.176 


0.388 


Ad 


-0.564 


1.000 


0.301 


-0.758 


0.032 


Ad 


-0.732 


1.000 


-0.170 


0.548 


0.190 


Am 


-0.827 


0.301 


1.000 


-0.385 


-0.257 


Am 


0.577 


-0.170 


1.000 


0.228 


0.524 


Ad 


0.346 


-0.758 


-0.385 


1.000 


-0.196 


Ad 


-0.176 


0.548 


0.228 


1.000 


0.450 


As 


0.075 


0.032 


-0.257 


-0.196 


1.000 


As 


0.388 


0.190 


0.524 


0.450 


1.000 



TABLE XVII: The statistical (pstat) and systematic (psyst) correlations of the second moments of the extracted helicity 
distributions in the measured x-range. 



Pstat 


A'-',/ 






A'-V7 


A(2)s p... 










A'^'s 


A(2)m 


1.000 


-0.803 


-0.522 


0.277 


0.048 A^^^M 


1.000 


-0.847 


0.430 


-0.399 
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